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Preliminary. — Professor Stokes's Report to the British Association, 1862, with Outline 

of the Method. 

Section I. — Review of Previous Experiments and Criticism of Fresnel's. 

In his report to the British Association in 1862, Professor Stokes says : " The exact- 
ness of the spheroidal form, assigned by Huyghens to the sheet of the wave surface 
within Iceland spar, does not seem to have been tested to the same degree of rigour 
as the ordinary refraction of the ordinary ray; for the methods applied by Wollaston 
(Phil. Trans., 1802, p. 381) and Malus (Mem. de l'lnstitut Sav. Etran., torn, ii., p. 303, 
1811) for observing the extraordinary refraction can hardly bear comparison for 
exactness with the method of prismatic refraction adopted for the ordinary ray; and 
observations on the absolute velocities of propagation in different directions within 
biaxal crystals are almost wholly wanting. 

" This has long been recognized as a desideratum, and it has been suggested to 
employ for the purpose the displacement of fringes of interference. 

" It seems to me that a slight modification of the ordinary method of prismatic 
refraction would be more convenient and exact. 

" Let the crystal to be examined be cut, unless natural faces or planes of cleavage 
answer the purpose, so as to have two planes inclined at an angle suitable for the 
measure of refractions; there being at least two natural faces or cleavage planes left 
undestroyed, so as to permit of an exact measure of the directions of any artificial 
faces. The prism thus formed having been mounted as usual and placed in any 
azimuth, let the angle of incidence or emergence (according as the prism remains fixed 
or turns with the telescope) be measured by observing the light reflected from the 
surface, and likewise the deviation for several standard-fixed lines in the spectrum; 
each observation furnishes us with an angle of incidence and the corresponding angle 
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of emergence, tlie angle of the prism being known. For if <j> be the angle of incidence, 
D the deviation, i angle of prism, xjj angle of emergence, 

D = cj)-\-\jj — i 

" But without making any supposition as to the law of double refraction, or assuming 
anything beyond the truth of Huyghens's principle, which, following directly from the 
superposition of small motions, lies at the base of the whole theory of undulations, 
we may at once deduce from the directions of incidence and emergence the direction 
and velocity of propagation of the wave within the crystal. For if a plane wave be 
incident on any plane refracting surface, it follows directly from Huyghens's principle 
that the refracted wave or waves will be plane, and if <£ be the angle of incidence, 
(j> the inclination of the refracted wave to the surface, V the velocity of propagation in 
air, v the wave velocity in the crystal 

sin <jf> sin $>' 

"T~~r~* 

" And if (j> \jj' be the inclination of either refracted wave to the faces of our crystal 

prism we have 

v sin <j)=Y sin <j/ .......... (l) 

v sini//=Vsin t// . , . (2) 

adding and subtracting (1) and (2) and remembering (3), we get respectively 

. + i/r 6— 4r TT . i <b f — i|r' ,. x 

V Sin ~ zr 1 - COS r — r J ~= V Sin ~ COS ^-TT V 4 ) 

2 2 2 2 v ' 



" By division 



or 



v cos —~z~^ sin ^rr- = V cos - sin x — ^- L ~ ...... (5) 

tan ?_— i- cot 2_jl == tan - cot — - — 

^ A Z & 

tan — ~^== tan - tan cot — ^ ......... (6) 

" Equations (3) and (6) determine <j> and t//, and v is known from 

sin </> ___ sin <£' 

or 

sin ^ sin ty' 

V v 
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" Hence we can find the velocity of propagation of the wave, the normal to which 
lies in a plane perpendicular to the faces of the prism, and makes known angles with 
those faces, and hence with the crystallographic axes." 

In accordance with these suggestions, I undertook a series of observations at the 
Cavendish Laboratory, Cambridge, which I propose in the present paper to describe ; I 
also wish to discuss the results arrived at, and to compare them with those deduced 
from Fresnei/s and some of the rival theories of double refraction. 

But previous to this it seems natural to devote some space to the consideration of 
the experiments that have already been made with the same object. 

These we may divide into two classes: (1.) Those which have reference to Iceland 
spar and other uniaxal crystals ; (2.) Those in which biaxal crystals were used, 

Brewster (Thirteenth Report of the British Association) proved conclusively that 
one wave in Iceland spar obeys strictly the ordinary law of refraction. 

Swan (Edinburgh Trans., vol. xvi., p. 375) obtained by direct measurements with 
prisms placed in the position of minimum deviation values of the ordinary refractive 
index, which differ at most by '00002. 

To verify his construction for the extraordinary ray, Huyghens himself made but 
few experiments, and it was not till 1802 that Wollaston (Phil. Trans., 1802, p. 381) 
attempted to test the theory with any degree of exactness. 

In 1810 Malus ("Theorie de la Double Refraction," Paris Mem. des Savants 
Strangers, torn, ii., p. 303) undertook a series of experiments with the same object in 
view. 

Rather later, Biot undertook the same task, while more recently Rudberg (Pogg. 
Ann., vol. xiv., p. 45) and Mascart have measured the values of the principal indices 
by means of prisms cut parallel to the axis. 

Since 1862 Professor Stokes has applied the method indicated above to prisms of 
Iceland spar, and finds the results of experiment agree with Htjyghens's construction 
with a possible difference of '0001 in the values of the refractive index. The details 
of the experiments are as yet unpublished. 

(2.) We must consider next the case of biaxal crystals. Various experimenters have 
determined the values of the principal refractive indices for different crystals. Fresnel 
alone has endeavoured to verify his theory by experiment ( f (Euvres Completes de 
Fresnel/ torn, ii., p. 415; second supplement to the " Memoire sur la Double Refrac- 
tion"). The method adopted was to observe the displacement of the fringes of inter- 
ference, formed from two parallel slits, produced by introducing two plates of topaz of 
the same thickness cut in different directions from the same crystal into the paths of 
the pencils proceeding from the two slits respectively, and to compare this with the 
displacement calculated on Fresnei/s theory. 

The first set of observations combined with the angle between the optic axes serves 
to determine the principal velocities a, b, c. a being assumed to be unity, let 

mdocclxxix. 2 p 
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C 2 =a?— y=l — y 

Then 

=0*00338 

y =0*01222 

y— /3=0*00884 

In the second set of observations the plates are inclined to the incident light, so that 
the light when in the crystal passes along an optic axis. 

To calculate the displacement of the fringes on Fresnei/s theory, it becomes neces- 
sary to know the velocity of the light in the crystal ; this is given by 

v^=l— -y sin 2 ^ 

X being the angle between the direction of vibration and the normal to that circular 
section of the wave surface whose radius is a or unity. 
Hence 

neglecting y 2 sin 4 x> & c - 

Theory gives for the displacement measured in wave lengths 

N=12-73 
Experiment gives 

N=11'87 

For the other two waves theory gives 

N=46'03 
Experiment gives 

N=45*49 

In the third experiment the light was incident at about 60°. 

In this case, tf> being the angle between the direction of vibration and the normal to 
the circular section radius b 

v*^— 1 ~/3 cos 2 <£-— y sin 2 <f> 
and to the same approximation 

i; l= i — -|(/3 cos 2 cj)-\-y sin 2 <f>) 

Whence theory gives 

N=19*57 
Experiment gives 

N=20'63 

Thus in each case the agreement between theory and experiment is fairly close ; but 
we must remember that the quantities to be measured are very small, and that in 
obtaining the theoretical results, small quantities have been neglected, without showing 
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the effects they would produce, and which, though small, may be sufficient to constitute 
the difference between Fresnel's and some rival theory. Moreover, the values of the 
principal refractive indices in topaz are so nearly equal that any conceivable wave- 
surface will differ but little from two spheres. 

To estimate then the weight which we may attach to these results, let us see how 
tar the experiments are consistent with one of the rival theories. 

Now the equation to the surface of wave slowness in the theories of Green and 
Cauchy involves more constants than the three principal refractive indices of the 
crystal. 

To determine these we must make the principal sections agree with experimental 
results in more points than the extremities of the axes ; and it is probable, therefore, 
that they will differ but little throughout. 

In fact, the existence of these unnecessary constants is a radical defect in both these 
theories. 

Lord Bayleigh, however, has proposed a theory in which the only constants are 
the principal velocities (Phil. Mag., vol. 41, series 4, 1871). 

The shape of the surface of wave slowness is also considerably different from that 
given by Fbesnel's theory. ... 

Its principal sections consist of a circle and the inverse of an ellipse instead of a 
circle and an ellipse. 

He supposes " the density of the ether in a crystal to be a function of the direction 
of displacement/' while the forces called into play are the same as in a homogeneous 
medium. 

From this he deduces as the equation to determine the velocity of a wave front, of 
which the direction cosines of the normal are Z, m, n 9 the equation 

aH 2 , 5 3 m 2 , cV 



V 2 — a 2 V 2 — 5 3 V 2 — $ 



a, 6, c being the principal velocities. 

Let us find the velocity for a wave perpendicular to the plane zx, making an angle 
a with Ox, l=- — sin a, m=0, n= cos a. 



V 2 (c 2 cos 2 a+a 2 sin 2 a)=a 2 c 



2 



or since a= 1, c 2 = 1 — y 



neglecting y 2 , &c. 



V 2 {(1 — y) cos 2 a+ sin 2 a} = l— ■ y 
V 2 = (1 — y)(l —y cos 2 a)"" 1 
V=(l— ^y)(l + iy cos 2 a) 

V= 1 — ^y(l — cos 2 a) 

= 1 — Jy sin 2 a 

2 o 
1 Zf 
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But replacing a by ^, this is exactly the form found above for v l in Fresnei/s 
second experiment. 

Now let us take a wave perpendicular to the plane yz inclined at an angle <fi to Oy. 
Then Z=0, m= sin </>, n= cos </>, and the equation for V 3 becomes 

5 3 sin 3 cf) $ cos 3 (p^ 

Y 2 —b^ V 3 — c 3 

V 2 (& 2 sin 2 <£+c 2 cos 3 <£)=& 2 c 2 
or 

V 2 (l~/3sin 2 (^--7COs 3 </))==(l~y8)(l--- r ) = l--( i 8+y) 
to the same approximation 

v= { i -W+y) } i i +W ^i>+r cos 3 <£) } 

= 1 -~J(/J cos 2 </>+y sin 2 </>) 

and this agrees with the formula used in Fresnei/s third experiment. 

Thus Fresnel's experiments afford as mtich a verification of Lord Rayleigh's 
theory as of his own, and are, therefore, an insufficient test of the truth of either. 

They are, however, the only attempts made to verify the theory. 

Section II. — Description of the Apparatus used and Method of making the Experiments. 

My own work was undertaken in the hope of obtaining results sufficiently accurate 
to decide between some of the various theories which have been propounded. 

I proceed to describe the apparatus used. 

The crystal itself was a piece of arragonite, very clear and of a light colour, obtained 
from Germany by the Demonstrator at the Cavendish Laboratory, Cambridge, through 
Hilger, of Tottenham Court Road, London, W.C. 

Originally it was in the form of an hexagonal prism. 

Four of the faces were those marked m in Professor Miller's mineralogy, the other 
two were the faces marked a. 

The natural ends had been cut off, and the artificial ends were approximately per- 
pendicular to the m and a faces. 

The crystal consisted of several twins, the twin planes being parallel to two of the 
m faces. 

At this stage Professor Stokes kindly undertook to examine the crystal, to deter- 
mine in which direction to cut it with the greatest advantage for the purposes of the 
experiment. He found that fair reflexions were obtained from three of the m faces, 
and decided that they should be left untouched to determine the position of any other 
cut faces. 

Again, the mean axis of elasticity bisects the angle between the contiguous m faces ; 
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the optic axes lie in the plane perpendicular to this axis. A prism formed with its 
refracting edge nearly parallel to the axis of b, would permit of observations being 
made in a plane passing nearly through the optic axes ; and it was found that, by 
inclining the faces so that the axis of a nearly bisected the angle of the prism, and 
making that angle about 42°, light could be passed through the crystal in the direction 
of either optic axis. It was therefore thought right that the prism should be so cut. 

Moreover, by polishing one of the a faces and cutting a face nearly coincident with 
the fourth m face, a second prism was formed with its refracting edge nearly parallel 
to c. 

It was important to get observations in different zones from the same piece of 
arragonite, for its chemical constitution is more or less variable, and the constants — i.e., 
the principal refractive indices — may vary in different specimens. 

The work of cutting was successfully performed by Hilger, and the faces of the 
prisms were very fairly plane. 

The crystal being thus cut, it remained to determine accurately the position of the 
cut faces with reference to the crystallographic axes. 

For this purpose, and throughout the observations, I used a goniometer, by Grtjbb, 
kindly lent me by Professor Stokes. 

The collimator is fixed to a graduated circle of about 9 inches diameter. The 
graduations, which are in silver, are on the flat face, not the limb of the circle, and the 
circle is divided to arcs of 10'. 

Attached to the reading telescope there are two verniers ; the graduations on the 
verniers being in the same plane as those of the circle ; and by means of these the 
circle can be read to angles of 10". 

In the centre of this circle there is a table to hold the prism; the table being 
attached to a second graduated circle with verniers which can be read to half minutes. 

It was requisite throughout the experiments to have the edge of the prism parallel 
to the axis of revolution of the telescope, 

To attain this the prism was attached to a stand, adjustable with set screws, which 
rested in slots on the table of the goniometer. 

The reading telescope was fitted with a needle point instead of the usual cross wires. 

In order to level the prism, I covered half the slit, so that when viewed directly the 
line of junction of the covered and uncovered portions coincided with the needle 
point. 

I then moved the telescope until the images of the slit, formed by reflexion at the 
faces of the prism, were brought in turn into the field ; and altered the levelling 
screws till the same coincidence was obtained. 

Homogeneous light was secured by using a Bunsen's burner with a platinum wick, 
one end of which was immersed in a saturated solution of salt-and- water. 

By the advice of Professor Stokes, the observations were made at angles of incidence 
increasing in arithmetic progression. 
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The direct reading was taken several times during each day s work to avoid error. 

The method of observing was as follows :— I set the telescope to a given scale 
reading, and turned the prism till the image of the slit seen by reflexion at one of its 
faces coincided with the needle point. 

The difference between the direct reading and the known one to which the telescope 
was set gave the supplement of twice the angle of incidence, whence the angle of 
incidence is easily obtained. I then turned the telescope till one image of the slit seen 
by refraction through the prism coincided with the needle point and took the reading. 

The difference between this and the direct reading gave the deviation for that wave. 
The observation was then repeated for the second wave. 

I then reset the telescope, so that its axis made an angle of 4° with it's former 
position, and took a similar set of observations. This gave a series of deviations 
corresponding to angles of incidence increasing by 2°. 

Just in the neighbourhood of the optic axes, observations were made at every 1° of 
incidence. 

For the first prism, each observation in the set finally chosen was repeated three 
times, and the mean of the results taken. 

The difference between this mean and any observation was only in a few cases as 
great as 10". 

In the case of the second prism, the observed values of the deviation on different 
days for the same angle of incidence agreed so closely, rarely differing by more than 
20", that it seemed unnecessary to take a third set. 

I found also that the refracted images were much more distinct for angles of 
incidence greater than that for minimum deviation, than for angles less. 

On passing through the position of minimum deviation, I therefore reversed the 
prism, so that the face of incidence became that of emergence. 

Fig. 1. 
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Thus, if A B C D be the ray undergoing minimum deviation, for a ray such as 
P B Q R, I made P B the incident ray, Q II the emergent ; while for a ray like 
S B T Uj 1 reversed the prism, so that U T became the incident, B S the emergent ray. 

The cause of the increased distinctness lies in the fact that as the angle of incidence 
increases from that giving minimum deviation, the breadth of the image of the slit 
decreases, and vice versd. 
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Section III. — Determination of the Position of the Principal Plane of the First Prism 
with reference to the Crystallographic Axes. — Values of the Quantities observed in 
the Experiments, and Calculation of the Reciprocal of the Wave Velocity for the 
First Prism. 

We proceed now to the measurements made to determine the position of the 
principal plane of the first prism, of which the faces are P P /# 

Fig. 2, 




Let A B C be the points in which the axes of the crystal cut a sphere, with its 
centre at the origin. 

Let M M' be the poles of the m faces of the crystal, M M' lie in the great circle 
A B, and the angle M M' is known; also B M=B M'; let P P, be the poles of the 
cut faces of the prism. 

Let 

PM=6>, PM'=0' 
P / M=^,PM'=r 

Let a, /3, y be the angles which P makes with the axes. 

Let MA=ju. 

From triangle PAM, 

cos PM= cos AP cos AM+ sin AP sin AM cos PAM 



or 



From P A M'. 



cos 0= cos ol cos /x+sin a sin /jl cos A 

cos 9'=z — cos a cos /x+ sin a sin /x cos A 
.\ cos #— cos a cos /x= cos $'+ cos a cos ft, 

cos 6 — cos 9 f 



* • OuS ot*— *» 



2 COS }M 



fO A JLJL j~ K-5 JL J. JL ,*<. 

COS fl v ' 
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Again, from triangle PMB, 

cos PM= cos BP cos BM+ sin BP sin BM cos PBM 
or 

cos 0= cos /5 sin /x+ sin f3 cos /jl cos B 
and from P B' M, 

cos 0'= cos /3 sin //,— sin /3 cos /x cos B 

• » COS o — • ~ ~ ~~ 

zsm/i ■ 

8 -\-0 8 — 8 
cos — — cos — — 

2 ... (2) 



* . • * . » 



sin /x, 



Now the angles 0' were determined by observation. 

Each angle was observed six times on two different days. 

The mean of the results was, 

6= 77° 18' 44" 

^=100° 43' 30" 

Some difficulty was^ experienced about the angle /*. 

Professor Miller gives it as 58° 5' ; but this, when substituted in the above formulae, 
gives a value for P V / the supplement of the angle of the crystal, which differs con- 
siderably from the measured angle. 

I therefore measured the angle /jl and found it to be 58° 30' nearly. 

Mr. Garnett (the Demonstrator at the Cavendish Laboratory) kindly communicated 
this result to the Rev. H. P. Gurney, of Clare College, who was then lecturing for 
Professor Miller, and he measured the angle with very nearly the same result as I 
had already arrived at. Taking, then, the value he obtained, and assuming as before 
that B bisects the angle M M', we have 

/x=58° 28V" 

Whence from the formulae (1) (2) 

ot=67° 10' 35" (3) 

(3=88° 52' 40" .......... (4) 

For the position of P, we have similarly, 



= 80° 39' 

0/=lOl° 16' 23" 
Whence, 

OL ZTIZQij t/ ^J\J 6 e « . » e . . . . \ *"W 

fj ^^ t/I DOU • e « • • . • • . • \ ^ / 

* This value of the angle ja is given in Dana's * Mineral ogy,* for specimens of arragonite coming from 
Silesia, The piece used -was obtained from that neighbourhood. 



MR. R. T. GLAZEBROOK ON PLANE WAVES IN A BIAXAL CRYSTAL. 

To determine the angles B A P, B A P /? the triangles BAP and B' A P y give 

(since AB=~ 

cos/3 
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whence 



similarly from B' A P y 



COS PAB: 



Sina 



BAP=88° 46' 40" 

B'AP=88°49' 10" 



Fig. 3. 




From triangle PAP, 

cos PP,= cos AP, cos AP+ sin AP, sin AP cos PAP, 

Whence 

PP=137° 9' 55" 

The mean of numerous experiments gave 

PP=137°9 / 30" 

The smallness of the difference affords a strong presumption in favour of the 
measurements. 

Let P P cut AC in L. 
Then 

sin Arr= . ^ sm AP 

' sm PP, 

Whence 

AP.P^O 3' 23" (7) 

Again 

cot LP / sin AP,= cos AP, cos AP L+ sin AP,L cot LAP, 

.\ LP =67° 27' 25". . (8) 

a -r t-> sin AP, sin LAP, 

sm ALP= t^-t- -' 

sin Lr / 

ALP=1° 12' 4" (9) 

MDCCCLXXIX. 2 Q 
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sin APsinAP.L 



sin AL : 



sin ALP 



/ 



AL=2° 31' 58" ' . . . . (10) 

The solution of the triangle ALP leads to the same results. 

We now come to the experimental determination of the deviation and the calcula- 
tions of the refractive indices. 

Throughout the work for this prism, <j> represents the angle between the normal to 
the wave incident on or emerging from the • face P, when in air, <f> the corresponding 
angle in the crystal, while \jj t// are the angles for the face P. I propose then to give a 
table of the values of (/> or \jj according as the wave was incident on P, or P. 

(2.) Of D+£, i being the angle of the prism, D+i is given because the expression 

, which is equal to — — , occurs in the formula, and it saved trouble to calculate 

D+i at once from the readings. 

(3.) Of the values of <// deduced from the formulae 

tan r , = tan - tan J —i~~ cot - 
2 2 2/ 2 

proved in the first section and ^'+^'=i. 

(4.) Of the values of fi the refractive index as the mean of the results, given by 

the two formulae 

sin </>=/x sine// 

sin \jj=z[jl sin \jj' 

Any difference in the two values of /x found from these two formulae was due only to 
the errors necessarily arising from the use of proportional parts, and rarely exceeded 

•00001. 

I also introduce a list of errors in the calculated quantities which would occur from 

errors of 10" in each of the observed quantities, taken so as to produce the maximum 
effect in the result. 

This list was obtained by setting down from the tables the errors to which the 
supposed errors in D and <f> would give rise. 

The numbers in the table of errors for jjl are the digits in the fifth place of decimals. 

The next page gives the complete work for one observation as a sample. 
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Work for reducing one Observation. 



Direct reading . 
Reflexion reading 
Deviation reading 
iz= angle of prism 



• • 



» • » 



• • 



* * • 



» • 



263° 39' 
180° 5' 
297° 34' 12 
42° 52' 30 



// 



// 



(j)= 1{ 1 80 — (Direct— Reflexion) } 
=4(180° 5'— 83° 39') 
=48° 13'. Error- 5". 



D+i=D-220° 48' 30" 

= 76° 45' 42". Error+10". 



xfj=J)+i-cl>=28 32' 42". Error+15". 
t±t=^±i== 38° 22' 51". Error+ 5'. 

A A 



^-^=^1^-^=9° 50' 9". Error-10". 

A A 



% 



L tan -=9*5936354 

A 



T x^+Y" 

Licot o 
10-1012497 


Error. 
216 


L tan 

9-2389842 


Error. 

-1252 


L tan ^ ^ 

8-9338693 


Error. 
-1036 


0'— -f' 

2 
4° 54' 30" 


—4" 


0' 

26° 19' 45" 


-4" 


16° 30' 45" 


+ 4" 


L sin 

9-8725466 


-94 


L sin 0' 

' 9-6469205 


-170 


log /a 

•2256261 


-76 


L sin ty 

9-6792905 


580 


L sin t^' 
9-4536616 


284 


logyll, 

•2256289 


396 


yii, from 
1-68123 


-3 


yii, from t^ 

1-68124 


+ 15 


mean 
1-681235 





Q A 
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Section IV. — -Description of the Second Prism: Position of its Axes. — Observed Values 
of the Deviations and Angles of Incidence, and Calculation of the Reciprocals of 
the Wave Velocity. 

It will be remembered that, in describing the cutting of the crystal, I stated that a 
second prism was formed with its edge nearly parallel to C, and its faces not differing 
much from the m and a faces of the crystal. 

Fig. 4 




Let the figure, fig. 4, represent a section of the crystal by a plane perpendicular to 
the m faces. 

The traces of the faces of the second prism are A B and C D. 

The greater part of the face A B lies in the twin crystal, so that the results for the 
second prism apply to that crystal. 

On examination, the m face B C was seen to consist of two parts, which are not in 
the same plane* 

Observation showed that the upper half (the prism being placed as in the figure) is 
in the same zone as the faces DE 3 E F, while the lower portion is considerably 
removed from this zone, 

I found also that when the faces D E, E F, were levelled, the line of junction of 
the dark and bright portions of the slit, when seen by reflexion from A B, coincided 
with the needle point ; so that A B is in the same zone as the m faces ; while the 
reflexion from D, though in the field, was too high. Thus if Q Q y , fig. 5, be the 
points where the normal to A B, fig. 4, and that to D, fig. 4 ; produced backwards 
meet the sphere of reference, A, B, 0, fig. 5, being the poles of the principal planes ; 
Q lies on the great circle AMB; Q / is slightly below it. 

Let M, M/, fig. 5, be the points where the normal to E F, and that to E D pro- 
duced backwards cut the sphere. Then 

AM=AM/=58°28'. 

The mean of several observations gave for MQ=103° 58\ Hence 

AQ=:45° 30' ........... (1) 

The mean of the observations for M Q / and M/ Q /? gave 
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MQ/=57° 51' 3O"=0say;. 
M/Q / =59°9 , = 0'say ; 



and the triangles MAQ / M/A Q, give, 

cos AQ/ 



cos 0-fcos 0' 

2 COS yL6 



where /x=AM=M/A. 

"Whence substituting for 6 and 0' 

AQ,= 2°43' 20" . 

Fig. 5. 



* • 



• *•»•• I ^ J 




In the triangle M A Q, the three sides are known ; we can thence find the angle 
MAQ / and obtain 

MAQ,= 76° 24" 45'. ......... (3) 

In the triangle Q A Q / the sides A Q, A Q / and the angle Q A Q / are known 

AQ = 45° 30' 
AQ /= = 2° 43' 20 
QAQ^IOS 35' 15 

Whence from the formula 



// 



// 



we get 



cos QQ,= cos AQ cos AQ/4- sin AQ sin AQ, cos QAQ / 

QQ / =46° 12' (or rather less). 



The mean of several experiments had already given for the value of QQ, 



O i r\' 



46° 10 

These results, differing as they do by less than 2\ afford a verification of the 
accuracy of the measurements. 
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The position of the plane Q Q / is completely determined by the arc A Q and the 
angle A Q Q // ; we have already found, 

AQ— 45° 30' 
For the angle A Q Q / the triangle Q A Q / gives 

sin AQQ,= ^qq' sin QAQ, 

whence 

AQQ,=4 o 37'20" (4) 

Thus we see that the principal section of the second prism is inclined to the plane 
c at an angle of 4° 37' 20". 

I mentioned above that the second prism lay in the twin crystal. This was proved 
by passing light through the crystal, in a plane nearly parallel to the principal plane 
B C, when four images of the slit appeared. On cutting off the end containing the 
second prism, two images vanished ; on cutting off the other end, only the second pair 
vanished, and the first pair reappeared. The work for the first prism was done with 
the second pair, the first pair being rarely sufficiently distinct for measurement. 

Throughout the work, <j> <j/ denote the angles which the normals to the incident or 
emergent wave, and the wave in the crystal make with the normal to the face Q. 

The same remarks as to the reversal of the crystal hold as in the case of the first 
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Section V. — Determination of the Values of the Principal Indices. — Angle between 

the Optic Axes. 

Our next step is to determine the values of the refractive indices on Fresnel's 
theory. 

It has been shown (Section III.) that the principal plane of the first prism cuts the 
plane C A C in a line inclined at an angle of 2° 31/ 58" to O A, and makes with that 
plane an angle of 1° 12' 4". 

Let us consider the intersection of this plane and the plane B C, B' O. 

Fig. 6. 




It will be a line, M M', say, inclined at 1° 12', or rather less to C C 

Take O A, O B, C, as axes of x, y, z, and let /v, /^, fi be the principal refractive 
indices, and therefore the semi- axes of the surface of wave slowness. 

The sections of the surface of wave slowness on Fresnei/s theory by the plane 
B C are a circle of radius \i a and an ellipse axes /^, [i in the directions of O C, O B 
respectively. 

Hence the values of the radius vector in this direction given by experiment — i.e., 
the values of the refractive indices for light traversing the crystal in this direction — 
are /x«, and the radius vector of an ellipse axes /x#, jju c , inclined at about 1° 12' to \L h . 

To determine which of the experimental values of /x is in this direction we require 
to find L M, L being the point of intersection of M M' and A C. 

From the triangle C L M, since C is a right angle, 

cos C'LM^ tan OX cot LM 
Whence 

LM=87°28' 4" ........... (1) 

And if P / be as before normal to the face P / of the first prism, 

LP =67° 27' 25" [Sebtion III. (8)]. 

Whence 

MP, =20° 0' 39" .......... (2) 
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and MP, measures the angle between the wave normal required in the crystal and the 
normal to the face P y . 

This is the angle we have denoted by <f>. 

From Tables II., III., VIII., and IX. we have corresponding values of 





4>' 


/* 




19° 32' 21" 


1-68134 




20° 34' 27" 


1-68130 




20° 28' 9" 


1-68119 




19° 29' 32" 


1-68123 


observations gave 








f 


fj. 




19° 39' 22" 


1-68126 




20° 41' 22" 


1-68129 



Again, none of the experimental values of /jl in Table III., lines 1 to 11, and 
Table IX., differ greatly from 1*68125, and this is not far from the mean of the six 
observations recorded above. 

We take, then, this as the radius vector of an ellipse axes /^, /x inclined at 1° 12' 

to /JL^ 

To determine /jl c we must have recourse to the second prism. 

Did the plane of the second prism coincide with that of x y, the values of /jl in 
Table XXI. would all be equal to /jl c . 

We shall consider later the effect of the inclination of the plane of the prism to that 
of x y, but for the present may remark that the values in Table XXI. differ little from 
1-53013, the value of (jl c found by Rudberg. 

Let us take therefore 

We have then in the plane BOC 

1 cos 2 # sin 3 # 

being the angle between /x and /^. 

Now 

^=1-68125 

^=1-53013 

6=1° 12' 

Whence 

^= 1*68132 (4) 

The value found by Budbehg is 1-68157. 
The difference is considerable, being "00025. 
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We may note, however, that his values of ^ deduced from two different pieces of 
crystal differed by '0004. 

Now [M a is the value of the refractive index for the first prism for the wave E when 



O r\ f c\c\ f f 



<£'=20° 39 
The values of/* at about 30' on either side of this are 



1-68578 
1-68574 



The value 



Table VI, 1-68589 
Table XII, 1-68569 

j*«=l-68580 . . . 



(5) 



seems to represent these fairly, 
consistent results. 
Rudberg gives 



Let us take it at present, and see how far it gives 



/x fl =l-68589 



To verify these results let us find the angle between the optic axes. 
If 28' be this angle in the crystal 



r 



tan S'= 




i 



1 


1 


2 


H>a 


1 


1 



V 



9. ..2 



l/v wj 



Whence 



8'= 9° 4' 58 



// 



(6) 



and if 28 be the angle in air as seen through a face normal to the least axis of the 
crystal 

sin 8= fib sin S' 
whence 

8= 15° 23' 30" (7) 



Kibchoff found 15° 27' as the result of experiments. 
The value given by Rudbekg's values of the indices is 

8=15° r 

which differs from Kirchoff's result by 20'. whereas the value given above as the 
result of mv experiments differs only by 3' 30". 
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Section VI. — Calculation of the Reciprocal of the Velocity of Wave Propagation on 
Fuesnei/s Theory for the First Prism and Comparison with Experiment. — The 
same for Lord Rayleigh's Theory. 

We are now in a position to calculate the reciprocals of the velocity of wave 
propagation for the first prism on the theory of Fresnel. 

Fig. 7. 




Let O 0', fig. 7, be the points in which the optic axes cut a unit sphere. 
L M the trace of the principal plane of the prism cutting A (7 A' in L, B C B' 
in M. 

Let N be any point in the plane, in which a wave normal meets the unit sphere. 
Let 

C , LM=x=l° 12' [Section III. (9)]. 

cos OL 
cos 0= cos LN cos 0L+ sin LN sin OL cos y= — — cos (LN— X) 

if tan X=tan OL cos ^. 

Similarly 

COS 0'=: r— COS (LM — X') 

COS X x ' 

if tan X'= tan OX cos x- 

Now 

AO=90°-9° 4' 58" [Section V. (6)]. 

— 80° 55' 2" 

AL=2° 31' 58" [Section III. (10)]. 

KJjLjZ^Z / o JjO 4 « . . . » , . . (ly 

whence 
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\=78° 22' 55" (3) 

X'=96° 33' 5" • • • (4) 

Also P / being the point in which the normal to the face P / meets the sphere 

LN=LP,+P,N 
=LP+f 
LP,=67° 27' 25" [Section III. (8)]. 
LN— X=LP,+<£'— X 
=<f/—lO° 55' 30" (5) 

LN-X'=LP+f-X' 

= f — 29° 5' 40" (6) 

Again, if ^ w 2 are the velocities of propagation in the direction 6 9', we know that 

v 1 a +t> 2 a =a a +c 2 '— (a 3 — c 1 ) cos cos 0' 
V\~ v a 3 =(a a — c 2 ) sin 0' sin 0' 

a c having here their usual meanings of the greatest and least principal velocities. 

2v 1 2 =a a + c a -(a a -c a ) cos (0+6') 
2V=a a +c 2 -(a 3 -c a ) cos {6-6') 

Prom these formulae we can find v x v z , and thence the values of ^ /x 2 for 



Taking 













1 






1 * /} , 

^2 




f^a~ 


= 1-68580 




fX e - 


= 1-53013 


c 2 = 


1 

aw 


= •351876 





_ l 


= •427117 



Q1 O 

O J. aJ 
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The Tables XXV. and XXVI. give the results of theory for the first prism. 

Let us take the outer sheet first. 

In lines 1-12, Table XXV., and lines 1-10, Table XXVI, the agreement between 
theory and experiment is fairly close, but these series of experiments overlap and are 
in the immediate neighbourhood of an axis of the section. 

From this point onwards, we see that the experimental curve lies entirely inside 
that given by Fresnel's theory, and the difference between the two increases as we 
approach the neighbourhood of the optic axes, line 16, Table XXV., 14 and 15 
Table XXVI., while for the rest of Table XXV. the differences diminish. 

The tables extend from about 8° on one side of the axis C, to nearly 16° on the 
other. This is obvious, for between the optic axes 0^0\ beyond them 0+0' is very 
approximately equal to twice NC, 

Now let us consider the inner sheet, as before in lines 1-14, Table XXV., lines 1—11, 
Table XXVI., the agreement is fairly close ; but these observations extend from about 
6° on one side, to 7° on the other of the maximum radius vector of the section. As 
in the case of the outer sheet, the differences here increase as we approach the neigh- 
bourhood of the optic axes, reaching a maximum of '00044, line 15, Table XXVI. 

After passing the optic axes, the differences decrease, and for the rest of Table XXV. 
average about *00024, and in this case the curve given by experiment lies outside that 
given by theory. 

On examining the differences given in Table XXV. for the outer sheet, I noticed that 
theory and experiment agree closely up to the immediate neighbourhood of the optic 
axis, line 16 ; while in the neighbourhood of the optic axis, Table XXVI., line 15, the 
increase in the differences is much more gradual. 

In fact, if we compare the experimental values of /x 2 in Table XXV., lines 11-14, 
and Table XXVI., 12-15, for about the same values of 0-\-ff we find considerable 
differences. 

This led me to calculate the results of other sets of experiments covering the same 
ground, which had been already made. 

The following table gives the results of one such set, taken about two months 
previously ; during the interval the goniometer had been taken partially to pieces, and 
afterwards refocussed and set again. 

It seems worth while to give the whole set of observations which cover nearly the 
same ground as Table II. , and form therefore a test of the degree of accuracy of the 
results. 

The measurements were made at angles of incidence, nearly, but not quite, coincident 
with those of Table II. 

In order to compare results, so far as the wave E (inner sheet) is concerned, a 
correction by interpolation is necessary therefore. 

The values of <//, it will be observed, differ by about 5' from those in Table II. 

In the table, the first column gives the values of </>' ; the second the values of ^ ; 
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the third those of jjl (corrected by Interpolation for the wave E) ; and the fourth the 
values of ju, from Table III. ; while the last column gives the differences between the 
two sets. 

Comparison of values of /x deduced from two independent series of observations at an 
interval of about two months. 

Table XXVII. — Wave O. 





0'. 


fi. 




o ( II 




1 


16 27 16 


1-68131 


2 


17 32 15 


1-68127 


3 


18 36 18 


1*68128 


4 


19 39 22 


■ • 1-68126 


5 


20 41 22 


1-68129 


6 


21 42 8 


1-68120 


7 


22 41 41 


1-68128 


8 


23 39 55 


1-68125 


9 


24 36 48 


1-68116 


10 


25 32 5 


1-68127 


11 


26 25 51 


1-68116 


12 


27 17 37 


1-68130 


13 


28 8 10 


1-68116 


14 


29 32 20 


1-68127 


15 


28 56 42 


1-68106 



Table XXVII. continued. — Wave E s 







Values of /t* 


Values of fx from 








reduced to the values 


Table III. corresponding 


Excess of the first series 


0'. 


JUL. 


of 0' in Table II. by 
interpolation. 


to the same values of 0' 
as in previous column. 


over the second. 


o ' // 

29 19 59 


1-68108 


1-68115 






29 43 %1 


1-68066 


1-68075 


1-68073 


-00002 


30 6 16 


1-68029 


1-68038 


1-68032 


•00006 


30 28 56 


1-67991 


1*68000 


1-67996 


-00004 


31 12 9 


1-67905 


1-67914 


1-67918 


- -00004 


31 53 13 


1-67828 


1-67837 


1-67841 


- '00004 


32 32 10 


1-67743 


1-67752 


1*67756 


- -00004 


33 8 28 


1-67668 


1*67678 


1-67674 


•00004 


33 42 29 


1-67590 


1-67600 


1-67591 


•00009 


34 13 53 


1-67508 


1-67518 


1-67519 


-•00001 


34 42 22 


1-67450 


1-67461 


1-67450 


•00011 


35 8 27 


1-67368 


1-67380 


1-67388 


-•00008 



Taking the observations for the wave O first, we observe that they do not differ 
greatly from those given for the same values of <j>' in Table II., except in the case of 
those which correspond to lines 12, 13, 14, of Table III., at which points the new values 

2 S 2i 
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of /x x are greater than those given in Table III., and agree better with the observations 
in the neighbourhood of the other optic axis. 

Another set of measurements in the same neighbourhood gave as values of fc 2 



1-68124 
1-68123 



1-68131 
1-68125 



Hence, on the whole, it seems likely that the results of the observations given in 
Table III., lines 12, 13, 14, are too small, and that an error has been made in the 
measurement at that point, owing to the indistinctness of the images which nearly 
overlapped, so that we may fairly replace Table II., lines 12, 13, 14, by Table XXVII., 
lines 12, 13, 14, and instead of the results of Table XXV. referring to these we shall 
have 



H x Theory. 


fi x Experiment. 


Difference. 


1-68114 

1-68104 
1-68096 


1-68130 
1-68116 

1-68127 


•00016 
•00012 
•00031 



This modification renders the results similar in the neighbourhood of the two optic 
axes. 

Let us now consider the results given in the latter part of the table for the wave E. 

The two sets of observations there recorded were made under entirely different 
circumstances. The telescope and collimator of the goniometer were removed in the 
interval between them, and both reset and refocussed ; and yet in but one case is the 
difference between the results of the two observations as great as *0001, while the 
average difference, irrespective of sign, is about '00005. 

This comparison forms a strong test of the accuracy of the experiments, and the close 
agreement of the results seems to show that the error in them is at any rate not greater 
than '00005. 

It would appear that more trustworthy results might have been obtained by com- 
bining the series of results registered iix this last table with those discussed throughout. 

But experiments covering the whole ground embraced in the series we have chosen 
had not been made at the same time previously. 

Moreover, we must recollect that each value of <f> and D+^ given in Tables I., IV., 
&c., is the mean of three for the first prism and of two for the second. 

And, again, this comparison assures us that the difference between the results of 
experiment and theory amounting, as it does in some cases, to '0005 in the value of jjl, 
is far too great to have arisen from experimental errors or errors of observation. 

It seems, then, that assuming the position of the plane of the prism, and the 
values of the principal axes to be accurately determined, we may assert that in a 
central section of the wave surface inclined at a small angle to the plane of the optic 



MR. R. T. GLAZEBROOK ON PLANE WAVES IN A BIAXAL CRYSTAL. 317 

axes, there is considerable difference between Fresnel's theory and experiment : that 
the differences between the two are most marked in the neighbourhood of the optic 
axes, and amount there to '0005 about. For the outer sheet of the surface (except 
just in the neighbourhood of the principal axes), the theoretical values of the radius 
vector are uniformly greater than the experimental, while for the inner sheet the 
reverse is the case. 

In fact, the curves of section of the surface, as given by experiment, approach more 
nearly to one another than those of the surface of wave slowness on Fresnel's theory. 

Fig. 8. • 




In the above figure, fig. 8, the dotted line gives the result of experiments, while the 
strong line gives the form of the section on Fbesnel's theory. 

To compare the results with those given by Lord Rayleigh's theory, we have the 
equation to the surface, a, 6, o being the principal velocities, 



p 






m 2 






# 






-1 


+ 






+ 






V 3 


V s 


•1 


V 8 


1 


o 

a" 






& 3 






c 3 





= 



1 

Put V= -, r is a radius vector of the surface of wave slowness 
r 

a?P hhw? cV 

jz^v+ 1 z jv + i - c ^ = 

Let us suppose as a first approximation that ra=0, i.e., that the plane of the prism 
is coincident with a principal plane of the crystal, 

aH % ( 1 — ch s ) -j~ chv 1 ( 1 — a%* 2 ) = 

r % c^(? = aW + 2 <m 2 

f — — -j_— - 

Also 1= cos 9* n= sin 9 

9 being the angle between r and the maximum radius vector 
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(1 1\ 

rj I O O ' OXXX K/ 



1 



(? 



2*84192 



1 1 



c M a 



2 



— '50062. 



From these and the known values of 6, fi or r can be calculated. The table gives 
the result. 

Comparison between Lord Rayleigh's Theory and Experiment. 



0. 


fi Theory. 


fi Experiment, 


Excess of Theory 
oyer Experiment. 


O i u 

1 31 10 

8 25 33 

11 6 40 

14 9 48 

15 48 28 


1-68570 
1-68261 
1-68028 
1-67688 
1-67476 


1-68568 
1-68190 
1-67918 
1-67519 

1-67274 


-00002 
•00071 
•00110 
•00169 
•00202 



The differences between theory and experiment are so marked that it seems un- 
necessary to calculate the values of fju for smaller intervals in the values of 6. 

And though some of the apparent difference may be due to the error made in 
assuming the principal plane of the prism to coincide with one of the crystal, that 
cannot account for the whole ; for we have seen that in Fkesnel's surface the error 
made by the same assumption appears only in the fourth place of decimals, in the value 
of the refractive index, while the differences between Lord Rayleigh's theory and 
experiment show themselves in the third place, and tend to increase with 6. 

Thus it seems clear that Lord Rayleigh's theory will not account for the phe- 
nomena of double refraction in arragonifce. This result agrees with that arrived at by 
Professor Stokes for Iceland spar. 



Section VII. — Effect of Varying Constants of Theory and Position of Plane of Prism 

with reference to Axes of Crystal, 

It remains now to discuss the effects of variations in the values of the axes a, b, c. 

The values of c and b are determined directly from observation, and in the neigh- 
bourhood of these axes the agreement is most close. Any alteration in them would 
affect especially that part of the section which lies between the optic axes, and for 
which the differences between theory and experiment are least ; let us consider, there- 
fore, a variation in the value of a, 
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Such a variation would produce an alteration in the angle between the optic axes, 
and therefore in the values of 0, ff for any wave normal, 

If, however, we take a wave normal at some distance from the optic axes, the values 
of 0, ff will be altered by nearly equal amounts in opposite directions, and 0+0' will 
be nearly constant. 

Let us then find the value for a in order that the experimental value in line 29 of 
Table XXV. may agree with theory. 

We have 

2^ 2 =« 2 +c 2 ~(a 2 --~c 2 ) cos (0+0') 

=a*{l-zo$(0+0 f )} 
+ c 2 {l+ cos (0+ 00} 

999 O + O 

Vf— <rcos^ — - — 

,2 fL 



.'. a" 



. «e+e f 

sm^ — — 

i 



v 1 =- 



/*! = l-67274 

0+0'=31° 45' 40 
c ? =-351876 



// 



From these we get 



a 2 ='447553 



as against '427117, and the value of fi c is 

^=1-49478 
instead of 1*53013 

A reference to the tables for the second prism shows at once the impossibility of 
such a change. 

It is true that any decrease in [x c would produce a decrease in the angle between 
the optic axes. This would increase 0' and decrease by nearly equal amounts, and 
on the whole produce a small decrease in 0+ 0'. 

This change would decrease v x and therefore increase fi l9 and hence to produce the 
required change in fx x — i.e., to make /x x equal to the experimental value in line 29 — we 
must take a value for \l c larger, but only slightly larger, than that given above. 

There remains now for discussion the effect of a variation in the position of the 
plane of the prism with reference to the crystallographic axes. 

We shall consider this by discussing separately the effect of a rotation about each 
of three lines through the centre approximately at right angles. These lines are 
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(1.) The normal to the principal plane of the prism ; (2.) The intersection of this 
plane with the plane B C B'; (3.) The intersection of the principal plane with the 
plane A C A'. 

(1.) Rotation about a normal to the principal plane of the prism. 

This cannot produce the required effect, for though we can thus decrease the 
differences on one side of the maximum radius vector of the section, we increase those 
on the other side by about the same amounts. 

(2.) Rotation about the intersection with B C B'. 

This is of no use, for the effects produced on opposite sides of the maximum radius 
vector are again of an opposite kind. 

(3.) There remains, therefore, only the rotation about the line of intersection with 
the plane C A C, or a variation in the angle x- This change is to be such that the 
sections of the outer and inner sheets may be brought into closer proximity; x must 
therefore be decreased. 

Now on referring to the tables of the experimental results, it will be seen that none 
of those in Table III., lines 1-11, Table VI., lines 16-29, and Table IX., lines 1-16, 
differ greatly from 1*6 8 125. 

If it were then permissible to neglect the angle ^, this would be the value of /x§, and 
the other section would be on Fresnel's theory an ellipse of fi a9 /jl c . 

In the abstract published in the 'Proceedings of the Royal Society/ No. 188, 1878, 
this has been done, and a limit assigned to the error thus introduced. 

Professor Stokes has since pointed out that near the optic axes the limit is con- 
siderably too small, and thus led me to the accurate calculation of the theoretical 
values of /x 1? /x 2 given above. 

The effect of a decrease in the angle x ^ay, therefore, be estimated by referring to 
the calculations for the case in which we put x= 0. 

The accompanying table gives the results. 6 being the angle between the wave 
normal and the major axis of the ellipse, which is given by 



0=f-0, 



o 



where Q is the angle between the normal to the face of the prism and major axip, 

and is 

^20° 0' 40" [Section V. (2)]. 



MR. R. T. GLAZBBROOK ON PLANE WAVES IN A BIAXAL CRYSTAL. 



321 



xaBTjE Jx.J\. V XJLJ.. 



Tables from 
which experi- 
mental result is 
taken. 


e. 


r or jx from Theory. 


r or jj, from 
Experiment. 


Excess of 
Experiment 
over Theory. 


jL.\. JL • 


8 25 38 


1-68193 


1-68204 




J? 


8 4 45 


1-68224 


1-68223 


- 1 


?? 


7 21 11 


1-68284 


1-68277 


- 7 


JJ 


6 35 31 


1-68342 


1-68333 


- 9 


JJ 


5 47 55 


1-68396 


1-68386 


-10 


?J 


4 58 26 


1-68444 


1-68439 


- 5 


5? 


4 7 22 


1-68486 


1-68480 


- 6 


V. 


3 42 55 


1-68504 


1-68500 


- 4 


XI. 


3 14 28 


1-68522 


1-68517 


- 5 


Y. 


2 38 14 


1-68542 


1-68540 


- 2 


XI. 


2 20 10 


1*68550 


1-68541 


- 9 


Y. 


1 34 26 


1-68566 


1-68570 


+ 4 


-r-ry 


1 23 23 


1-68570 


1-68560 


-10 


Y. 


0. 31 31 


1-68578 


. 1-68578 


- 


XI. 


27 12 


1-68578 


1-68569 


- 9 


Y. 


30 15 


1-68578 


1*68589 


+ 11 


-"cr-r 
-A.J.. 


31 18 


1-68577 


1-68574 


- 3 


y\.j.» 


1 31 3 


1-68567 


1-68570 


+ 3 


Y. 


1 31 10 


1-68567 


1-68568 


*x* JL ' 


j» 


2 30 43 


1-68544 


1-68560 


• • 


XL 


2 31 48 


• 1-68544 


1*68539 


- 5 


V. 


3 29 17 


1-68515 


1-68517 


+ 2 


XL 


3 33 40 


1-68510 


1-68502 


- 8 


V. 


4 26 32 


1-68472 


1-68467 


- 5 


J? 


5 22 18 


1-68422 


1-68421 


- 1 


J> 


6 16 31 


1-68367 


1-68376 


+ 9 


?} 


7 9 25 


1-68300 


1-68302 


+ 2 


J? 


8 38 


1-68230 


1-68226 


- 4 


?J 


8 25 33 


1-68194 


1-68190 


- 4 


JJ 


8 50 1 


1-68156 


1-68153 


- 3 


JL JL. 


9 14 14 


1-68118 


1-68131 


• • 


j? 


9 37 36 


1-68077 


1-68073 


- 4 


55 


10 40 


1-68037 


1-68032 


- 5 


j j 


10 23 7 


1-67997 


1-67996 


- 1 


55 


11 6 40 


1-67914 


1-67918 


+ 4 


55 


11 47 59 


1-67831 


1 '67841 


+ 10 


55 


12 27 8 


1-67748 


1-67756 


+ 8 


i i 


13 3 53 


V67666 


1-67674 


+ o 


55 


13 38 11 


1-67587 


1-67591 


+ 4 


?J 


14 9 48 


1-67511 


1-67519 


+ 8 


55 


14 38 44 


1-67439 


1-67450 


+ 11 


55 


15 4 51 


1-67373 


1-67388 


+ 15 


55 


15 28 8 


1-67310 


1-67328 


+ 18 




15 48 28 


1-67250 


1*67274 


+ 24 



MDCCCLXXIX. 






322 ME. E. T. GLAZEBEOOK ON" PLANE WAYES IN A BIAXAL OEYSTAL. 

For about 10° on either side of the major axis (the passage across which is indicated 
by the dark line in the table) the results of theory and experiment agree fairly well. 

The average difference irrespective of sign is (if we neglect the values for 6=2° 30' 43" 
and for $=z9° 14/ 14", which gives a point on the circle of radius 1*68125) '000047, 
and only in two cases does the difference amount to '0001. 

But when we come to values for 9 greater than 11° we see that the experimental 
value of fi is, as before, always greater than the theoretical, and that the difference 
increases with 9 and at last becomes nearly the same as in Table XXV. 

Thus this change does not produce agreement between the two. 

We may, however, consider shortly what alteration in /jl c would in this case bring 
theory and experiment more closely into agreement. 

We have 

1 ___ cos 3 sin 3 6 

"~"" " MM nJ llllil ' ' ' 

/^ fta, [Ac 

. \ Bjjl = — - smWSfic 

Taking the last values on Table XXVI II. 

S/jl== '00024 

0^15° 48' 

jul=1'67274 

/V-J>53013 

Whence 

S/^=-00248 

so that the minor axis of the ellipse, which having the same major axis as the experi- 
mental curve would pass through the extreme experimental point, is 

/V=l"53261 

instead of 1*53013 

And on referring to Table XXI., which gives the values of /jl very approximately, we 
see that this value is quite out of the question. 

Thus a decrease in the angle x will n °t produce the required effect. 

Neither will an increase. For the first part of the inner sheet in Table III., and 
throughout Table IX., we must have the theoretical value of jjl 1 and therefore of 6-^ff 
nearly constant; and this is clearly impossible along a plane section inclined at a finite 
angle to the plane of the optic axes. 
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Thus a change in ^ will not render theory and experiment consistent. 
Thus no change in the position of the plane of the prism will bring Fkesjstel's theory 
into agreement with experiment. 



Section VIII. — The Theoretical Investigations for the Second Prism on Fkesnei/s 

Theory with Comparison with Experiment. 

Our next step is the theoretical investigation for the second prism. 
Now experiment has shown that this prism lies in the twin crystal ; we must 
therefore find the axes of this crystal. 

Fig. 9. 




Let the figure represent a section of the crystal by a plane perpendicular to the axis 

oa 

The twinning takes place about an axis perpendicular to m\ and consists of a 
rotation through 180° about that axis. 

Its effect, therefore, is to bring B 5 fig. 9, into the position B 7 , where 

angle mOB = mOB' 

that is, to turn the axes in the direction B to A, through an angle, 

= 2(mOB)=2^-AOm) 

= ir~-2/x=180 o -~116° 56' 
= 63° 4' 

and if Q is the point where the perpendicular on the face Q of the second prism meets 
the sphere, we know that 

AQ=45° 30' [Section IV. (1)]. 

Hence since AA^GS 4' and Q lies in plane AOB, 

Now let A' B' C, fig. 10, be the points in which the new axes cut the sphere, and let 

2 T 2 
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PQM be the plane of the prism ; P being any point in the plane ; Q the normal to 
the face, and M the point in which the p]ane cuts B' C, 0, / the optic axes. 
Then if P is the pole of any wave in the crystal, PQ—^'. 

Fig. 10. 




Let PO=0, PO^^, we shall require the values of $, 6\ 
Now in triangle A 7 M Q 



angle Q=4° 37' 20" [Section IV. (4)]. 

A , Q=17° 34' 
cot MQ= cos Q cot A'Q 



whence 



MQ=17°3r 14" . . . . 



- (2) 



Let A' M Q=x- 

Then from triangle A' M Q. 

whence 



cos M= cos QA' sin Q. 



One/ Ar\'f 



x =85" 35' 40 



• » 



• * * 



. , \6) 



Now, 



put 



cos ff= cos PM cos MO'+ sin PM sin MO 7 cos x 
= cos M.O'(cos PM + sin PM tan MO' cos X ) 

tan MO' cos ^= tan X' 
ff _ cos MO' cos (PM— V) 

COo U •— — ■ ■ r~r • » » • » 

cos\ 



• • • ( 4 ) 
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Similarly 

* cos MO cos (PM + \) , . 

COS v — — ~ " . • * , , , , « , (0) 

cos X v f 

where 

tan\=tanMOcosx 

- Now 

^0=^0=90°- 9° 4' 58" [Section V. (6); 

= 80° 55' 2" 



From triangle A'MQ 



Whence 



sin A'M= sin MQ sin Q 



XL 1VI L JiO Z . , . (Q\ 

.-. MO=79°31' 10" (7) 

MO'=82° 18' 54" (8) 



Hence 



X — 22 oo 2 (9) 

X'=29° 39' 5" (10) 

PM+X=f +QM+X 

= f+40°]0' 16" (11) 

PM-X'=f+QM-X' 

Substituting these values of X, MO, PM+X, &c. (the values of tf being taken from 
Tables XIV., XVII., and XX.), in equations (4), (5), we can get a series of values 
for 6 ff. 

And as before [Section VI], remembering that & here is ir—ff in Section VI., 

2v 1 s =a*+c s +(a s -c s ) cos (d—ff) 
2V=^+c 3 +(a 3 -c 3 ) cos (6+6') 

Thus we obtain a series of values for v l9 % and hence for /x 1? /x 2 , given in 

m 1 i "V"VT"V 

JL cL L/J.t5 -ZV-x\- JL J . V » 

To obtain the distance of P from the minimum radius vector, we have to add to <£' 
the angle M Q=17° 37' 14" [Section VIII. (2)]. 
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Considering first the inner sheet, we see that theory and experiment agree closely. 

The greatest difference is '00013 ; and in this case, by referring to the adjacent 
experimental values, it is clear that the value 1*53003 is too small. 

We must remember, however, that the investigations are concerned with that sheet 
of the surface of wave slowness of which the section by the plane, x y, is a circle, and 
that the section considered is inclined at an angle less than 5° to this plane ; so that 
the form of the curve of section differs but little from a circle, and we may therefore 
reasonably suppose that any theory in which, one of the principal sections of the 
Surface of wave slowness is a circle, would also agree fairly closely with our results. 

We therefore proceed to discuss the results for the outer sheet given in the table 
just preceding. 

And here we see at once that there is considerable difference between theory and 
experiment. 

Moreover, from f =12° W 53" to f = 20° 5' 38", the theoretical values of /x are 
uniformly greater than the experimental; while from <£'=20° 35' 20" to cj)'=3l° 42' 34" 
the reverse is the case. 

The difference between the two increases on the whole uniformly, and in the last 
case is as great as "0005. 

For the values <f>= 19° 13' 4" and </> / =20° 5' 38", theory and experiment agree closely. 

To find the angle between any radius vector and the minimum one of the curve, we 
have to add to the value of </>' the angle M Q, which is 17° 37' 14". 

Let us call the angle Q O P, 9. 

Then we see that for the angles 0=0, 0=37° (about), and 0=90°, the curves given 
by theory and experiment agree. 

From 0=29° 54' 7" to 0=37° (about), the theoretical curve lies outside the 
experimental. 

From #=37° (about) to 0=49° 19' 48" the reverse is the case. 

In the neighbourhood of the axes the two must agree more closely. 

These conditions would all be satisfied by a curve which agrees with Fresnei/s 
section at the extremity of the minimum radius vector; lies inside it for about 40°; 
then cuts Fresnex/s curve again, and lies outside for the rest of the quadrant, agreeing 
again at the extremity of the maximum radius vector. 

At first sight there seems some discrepancy between this result and that arrived at 
already in Table XXV. for the first prism. 

But closer examination shows that the two results are corroborative. 

For, in the first place, the observations there tabulated apply to the neighbourhood 
of the major axis of the elliptic section. 

For about 6° on either side of this principal axis, the differences between theory and 
experiment are sometimes positive, sometimes negative ; that is, between these limits 
Fresnel's section satisfies the results of experiment; but in the last eighteen 
observations there recorded, the experimental value of /Xj is always greater than the 
theoretical. 
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Thus the experimental results are represented by a curve, which near the major 
axis agrees with Fresnei/s section, but which, as we go from 6° to 16° away from 
that axis, differs from that section, and lies outside the section. 

This result, viz. : that on going away from the maximum radius vector the curve 
given by experiment should lie outside that given by theory, is exactly that to which 
we have been led by the work for the second prism. 

The cutting of the crystal would not allow of observations being taken near the 
minor axis of the ellipse, in the case of the first prism. 

Of course, it has here been assumed that the sections of the surface of wave 
slowness, by the principal planes, are curves of the same kind ; or rather the results of 
experiment lead to that as a fact. 

Lord Rayleigh's theory needs no further discussion here, for his curve being the 
inverse of an ellipse, always lies outside the ellipse with the same axes, and therefore 
cannot agree with experiment. 

Section IX. — Possibility of an Error in the Position of the Crystallographic Axes 
discussed. — No Section of the Surface of Wave Slowness in Fresnei/s Theory can 
agree with the experimental Results. 

But the possibility remains that there is an error in the determination of the posi- 
tion of the faces of the prisms, especially in the second case, with respect to the 
crystallographic axes. 

The probability of any error is small, as the observations were repeated on several 
different occasions with nearly coincident results. 

To make certain, however, let us solve the inverse problem. Having given the two 
values, /x x /x 3 , of the radii vectores of the surface of wave slowness drawn in the same 
direction Imnm the crystal to find Z, m, n. 

We have as the equation to find v on Fresnei/s theory 

72 mo2 q$ 

L i JIl — i—^—o 

\ ..9, 12 T „o o v 



v*—a? v^—ffi v^—c 1 
Whence 

^_t^{P(&s+c s )+m s (c a +^ 

.\ if v l v % are the roots of this equation, that is, the velocities for the same wave front 

v^v 2 2 = PWtf-^mh^a 2 -^ n % a?W 
v *+ v *=z Z 3 (6 3 +c 2 ) +m 2 (c 3 + a 2 ) + w s (a*+ 6 s ) 

but n 3 = 1 — m 3 — P 
Whence we get 

&ap( c a _ a?) + a*m\£ - W) = v 1 V - <**&* 

P(c 2 -~-a 2 ) + m 2 (c 3 -~& 3 ) = ^ 3 +t' 2 3 -a 3 ~6 s 
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Whence 

Z 2 (a a — c 2 ) (a 2 — ¥) = (a 2 — vf 1 ) (a 2 — v^) 

m 2 (6 2 — c 2 ) (6 2 — a 2 ) = (6 2 — v^) (6 2 — ^ 3 2 ) 

n 2 (c 2 — a 2 ) (c 2 — 6 2 ) = (c 2 — ^ 1 2 ) (c 2 — Vc?) 

From the tables for the second prism, we see that for the inner sheet for 



For the outer sheet, 



<£'=28° 11' 21"- = 1"53021 



<b'=28° V 21"- = 1'68391 

r v 



Since the value of v l changes very slowly, we may reasonably treat these two values 
for <f> as coincident, and substitute in the above formulae. 
"We find 

cos" 1 ^=49° 35' 

cos"" 1 m 1 =40° 27' 
Again, for inner sheet for 



For outer sheet 



<4'=38° 2' 41" - = 1-53037 

v 1 



<£'=31 9 42' 34" - = 1*68437 



V 2 



Treating these two as coincident, we find 



For inner sheet for 



For outer sheet 



Whence 



cos"" 1 n 2 =55° 44' 
cos 1 " 1 m a = 3 4° 23' 



<i'=26 5' 11" - = 1-53028 



6'=26° 19' 10" -=1-68378 



cos- 1 n 3 =47° 56' 
cos- 1 w 3 =42° 8' 



We may justify the treating of the two values of <j> as the same, by the considera- 
tion that /x x varies at the rate of 00001 (about) for 1°, and therefore the difference for 
20' is inappreciable. 

Now the angle between the radii vectores in the first and second case cannot be 
less than 

MDCCCLXXIX. 2 TJ 



Q 
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cos" -1 n % — cos" 1 n x 

t'. c, 

not less than 55° 44 / -49° 35' 

or 6° 9' 

But the angle between these radii is the difference between the two corresponding 
values of <f>', or 

=^31° 42' 34"-28° V 21" 

= 3° 41' 13" 



Again, taking the second and third, the angle deduced from theory is not less than 

cos"" 1 %— cos"" 1 n 3 
or 55° 44 / — 47° 56' 
or 7° 48' 

31° 42' 34 // -26° 19' 10" 

= 5° 23' 24" 



but from experiment it is 



In both these cases the angle, as given by theory, must be very much larger than 
that given by experiment. 

Thus, in order that the three pairs of values of /jl 1 /jl 2 may correspond to three radii 
of the surface of wave slowness on Fbesnei/s theory, the angles between them must 
be half as large again as they are found to be. 

Hence no section of the surface of wave slowness on Fresnei/s theory can agree 
with the results of the experiments. We may remark, in addition, that the differences 
between consecutive values of /x for the outer sheet given by experiment are about 
double those given by theory. 

After this investigation it seems needless to discuss possible alterations in the posi- 
tion of the plane of the second prism, with a view to bringing the results of experi- 
ment more nearly into accordance with theory. While as regards the constants a, 6, c, 
it has already been seen that they have received their most probable values. 

Thus our results, so far as they go, point to the fact that Fresnei/s construction 
does not represent Nature, and that some other theory must be sought to explain the 
phenomena of double refraction. 

JL 

And as regards the rival theories already proposed, Lord Rayleigh's, we have seen, 
is even more than Fresnei/s at variance with observed appearances, while I can hardly 
think that sufficient data have yet been collected to repay the labour of comparison 
with the theories of Green and Cauchy, for on account of the undetermined constants 
it would be necessary to force theory and experiment to agree in so many points that 
they must, perforce, almost agree exactly along the small arcs investigated. 
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PART II. 

Section I* — Description of a Second Crystal and Results of Experiments. 

Before asking permission to lay the results obtained in the previous part of the 
paper before the Society, I thought it would be better, if possible, to repeat the inves- 
tigation with a second piece of the crystal. With some trouble I obtained a suitable 
piece last February from A. Hilger, of Tottenham Court Road. 

Like the former piece, it was an hexagonal prism in shape, but of considerably 

smaller cross-section. 

Fig. 1. 
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The axis of the prism was parallel to the axis of c, the base being perpendicular to 
the axis, while the top was broken off obliquely. 

The planes m / m gave the best reflexions, and were therefore chosen to determine 
the position of artificial faces. 

My chief aim in having the prism cut was to get as long a continuous arc on the 
wave surface as possible. To attain this I had a face R, fig. 2, polished as nearly 
coincident with c as possible. 

Pig. 2. 




The oblique end was cut (P, fig. 2) so as to be inclined at about 37° to R, the edge 
of the prism thus formed between P and R, being nearly parallel to the intersection of 
m and c. 

# N.B. The two Parts being quite distinct, the numbering of the Sections and Figures is commenced 
afresh with Part II. 

2 u 2 



332 MR. R. T. GLAZEBROOK ON PLANE WAVES IN A BIAXAL CRYSTAL. 

The face m / was cut so as to be inclined at about 35° 20' to P, the line of inter 
section of this face (Q, fig. 2) and P being also nearly parallel to that of P and R. 

So that the principal planes of the prisms formed by P and R, P and Q respectively, 
were nearly coincident, and as my observations extended almost from perpendicular 
incidence on Q to perpendicular incidence on R : they embraced an arc of over 70°. 

The instrument used and the methods adopted for levelling, reading, &c., have been 
fully described in the previous part of the paper. 

The results of the experiments are contained in Tables L, II. , III., and IV. 

Tables I. and II. refer to the prism R P. 

(/> is the angle made by the wave normal in air with the normal to P. 

if> the corresponding angle in the crystal. 

D is the observed deviation. 

i the angle of prism. 

D+i is given, for it occurred in the calculation, and was just as easily found as D. 

The value taken for i is 

i=sr 2' 56" 

This value is the mean of twelve observations, none of which differed from the mean 

by 10". 

<f> is found from the formulae <j&'+i//=& 

and tan - = tan - cot — - - tan ^~~ — 

which has already been proved (Part I., Section I.). 

In Table I. <f> is the same for both outer and inner sheet, being the angle of 
incidence. 

In Table II. <f> has different values for the two sheets, being the angle of emergence, 
for as already explained, on passing through the position of minimum deviation, the 
prism was reversed so that the face of incidence became that of emergence. 

Throughout the work //, x refers to the inner, /% to the outer sheet. 

The values of 8<fS, S/x are calculated as before on the assumption of errors of 10" in 
2c/> and D+i, the observed quantities taken; these errors being combined so as to 
produce the greatest possible effect in the result. 

Tables III. and IV. refer to P Q, for which ^=35° 18' 50". 

The symbols have the same meaning, but while in I. and II. the wave normal falls 
on the same side of the normal to P as the axis C, in III. and IV. it is on the 
other side of P. 

In III. <f> is the angle of emergence, and therefore different for the inner and outer 
sheets ; in IV. it is the angle of incidence, and the same for both. 

Each observation was repeated twice on different occasions; the results of two 
measurements rarely differed by 20'. 
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Section II. — Determination of the Position of the Planes of the Prism P IT with 

reference to the Axes of Elasticity. 

In order to compare these results with theory we require to know accurately the 
position of the principal planes of the two prisms. We will consider the prism, faces 
P It, first. 

Let A, B, O C be the directions of the principal axes of the wave surface, B 
bisects the angle between the? faces m m', C is parallel to the intersection of these 
faces. 

Let the normals to the faces m m' cut a unit sphere, centre O, in M M', the normal 
to the face P in P, and the normal to II in It. 

Let the direction angles of P and II with reference to A, B, C be a ± /J x y i9 a 2 /3 2 y 2 
respectively. 

[P lies in quadrant A B C, II in A 7 B C] 

Fig. 3. 




Let 

PM=^, PM'^0/ 

BM=0 a ,BM'=0 a ' 

Observations on the angles between P and m 9 P and m' give l9 0{ respectively. 
While observations on the angles between It and m, It and m give 0%, 0/ respec- 
tively. 

Let 

MM / =2^=2MB. 

Then, as in the first part, from triangles M' B P, M P B we get 

0/ -\- 0i 0i'--—0-i 
cos -^-z — - cos -^-z — - 

008/3] = . , (1) 

COS/4 x 7 



with similar formulae for /3 2 , a 2 . 

MDCCCLXXIX. 2 x 



\'^T0\ • 01 — 01 

sm -^-z — - sm -^-r — - 
cosa.= ~ (2) 
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All the angles Q x , 6{, Q % , 6 Z ', and yw were observed several times throughout the 
work. 

The mean of the observations give 

ft, =54° 17' 6" (3) 



ft/=75° 16' 20" . . 



ft =88° 40' 35 

ft-,' =88° 20' 20 
H =31° 53' 30 



// 



// 



// 



. . . (4) 
• • ■ (5) 

(6) 

(?) 



• * 



• ft 



In no case was the difference between the mean and any one of the observations so 
great as 20". 

Substituting these values we find 

a 1= 71 49' 30' 
ft=60 o 25' 49' 



,// 



J' 



. . . (10) 



xxlSO 



whence 



a 3 =90° 19' 10" 
/3. 2 =88° 14' 32" 

sin 3 ji= cos 3 aj+ cos 3 /S x 
sin 2 y ■%= cos 2 a 2 + cos 3 /3 3 

7l =35° 42' 57" . . . 



(11) 



77- r3 =l 47' 11 



// 



(12) 
(13) 



Let R O meet the sphere again in R', and let the figure represent the sphere as seen 
by an eye on O C produced. 

Fig. 4. 
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We have now sufficient data to calculate the angle R' P. 

CR':=7r-- • y 2 

AR'=7r— a s =89° 40' 50 

B>T>' . Q 
XV — P*9* 
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// 



From triangle B' A R' 



whence 



jd rom xriaiiofie o xL jt 



cos B'R'= cos BAR' sin AR' 

sin AR 

olli. Ota 



B'AR'= 88° 14' 32" 

.\ R , AC= 1° 45 ; 28 ;/ 

cos CP= cos CAP sin AP 



.\ cos CAP 



cosCP 
sin AP 



, *» OjAJl — oL 17 ^0 
.-. R'AP=33° 2' 48" 



From triangle R' A P 



cos RT= cos R'A cos AP+ sin R'A sin AP cos R'AP 
= cos (fir— a 3 ) cos a x + sin (77—- a 2 ) sin a x cos R'AP 



whence 



RT=37° 2' 53 



But RT can be measured directly, the mean of twelve observations, none of which 
differed from the mean by 10", gave 

RT=37° 2' 56" 

The close agreement between these results confirms the accuracy of the values of 
a l9 &, &c. 

Let P R / cut C A in L, we proceed to determine C L, R' L, and the angle ALP. 
From triangle AR'L 

cot AL sin AR'= cot AR'P sin R'AC+ cos AR' cos RAO 



.*. cot AL sin a 



^ cot ART sin R'AC+ cos a 2 cos RAG 

2 x 2 
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we have seen already that 



For A R' P we have 



R'AC=1°45'28 



// 



On substituting we get 



S1H -clJlv Jr — — ~ — — - gjj^ j^, ^fxx 

sm E P 
.\AB/P=59° 18' 37" 



AL — 88 08 18 

• » \J JLi — — 1 ^j L tm • * . 



* « * * • 1 10 1 



FT~k /Til i * 1 TPk/ AT" * 

or R L the triangle HAL gives 



sin sa> j. j — — 



sin R'AL sin AL 
sin AR'L 



whence 



R / L=2°2 / 37 



For ALP the triangle ALP gives 



. A T p sinAP . T A p 

sm LP 

LJr =Jl\/ x — Jt(jL=o5 19 



Whence 



• * » • » • • ■ I lOI 



ALP =5 9° 20' 11 



// 



♦ « ♦ * » 



• . . (17) 



Thus we have determined the position of the plane R' P completely. 



Section III. — Determination of the Position of the Plane P Q with reference 

to the Principal Planes. 

Our next step will be to determine the position of P Q% Q' being the point in which 
Q O produced backwards cuts the sphere. 

Fig. 5. 




Let the direction angles of Q' be a s /3 3 y 3 



Let 
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WQ'=0 a ' 



then 



cos /3c 



0/ + ft? 6 a 0<i 

COS 6 - 3 COS "^"77- 
COS /X 



tfq H~ 0O . Uo — 0O 

sin -^r-^ sm - * 



cos a 3 = 



COS fl 



But # 3 , # 3 ' are the angles between the faces Q and m, Q and m' respectively. 
And observation shows that 



3 =19 O 0' 15" 
0: 64°37 / 0" 



(3) 
• 3 b4 3/ u j 

Whence 

j3o — — O D Do 4i . « . » « • » . . . x^-j 

a 3 =60 o 42' 56" (5) 

JiA= - — ~r 

Q'BA=33° 37' 51" 
cos CQ'= sin BQ' sin Q'BA 

.•. y8 =CQ'=71° 0' 44" (6) 

To corroborate these values let us find P Q' and compare with experiment 

cos Q'AB= ^ 

sm a 3 



\Q'AB=21° 54' 12 



// 



But we know from the previous work (Section II.) that 

PAB=90°-PAC=58° 42' 42" 
.\PAQ'=36°48' 30" 

cos PQ'= cos AP cos AQ'+ sin AP sin AQ' cos PAQ 

Whence 

PQ'=35° 18' 35" 

The mean of a large series of experiments gave 

PQ'=35° 18' 50" 
and again we have strong evidence in favour of the correctness of the work. 
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Let the figure, as before, represent the sphere as seen by an eye on C produced. 

Fig. 6. 




Let Q' P produced cut C A in L' we proceed to determine C L' ; P L' ? and P L' A. 
From the triangle L' A P 

cot AL' sin AP= cot APL' sin L'AP+ cos AP' cos L'AP 

From triangle P Q' A. 

sin APL'= V-^x7 sin PAQ' 

sm PQ 

a— APL'=64° 41' 28" 
Also 

L'AP=CAP=31° XT 20" 

[Section II. from, previous work.] 

Substituting these values we get 



// 



The triangle AL'P gives 



AL / =88°44 / 20' 

CL'=1 15' 40" (7) 



sin ALT= . — sin APL' 

sm AL 



whence 

ALT=59° 13' 2" ......... . (8) 

Also 

sin L'P= - — T^zzr, sin L'AP 
sm APL 

whence 

To corroborate these results draw L K perpendicular to P L', then L K L' is 

approximately a right-angled triangle, and angle LL / K=60° nearly. 

Hence 

L'K = ^LL 7 approximately 

Now L'K^PL'— -PL^S' approximately 

LL'=CL — CL'^6' approximately 

.\L'K is equal to ^LL' 
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Section IV. — General Theory. 

To apply these results to the theoretical calculation. 
Taking the same figure 

Fig. 7. 




Let the optic axes cut the sphere in 0'; consider a plane cutting the plane A C 
between and 0' in L say; let N be any point on it; let NLA=x, NO=0, NO' = 0'. 

Then we know that if v l9 v% are the velocities along the wave normal through N, 
a, b, c the principal velocities, a being the greatest 



^ 1 S — ^ 2 2 =(a s — c s ) sin 9 sin 9' 
Vi*+v^=(a 2 +c 2 )--(a?---c 2 ) cos 9 cos 6' 
2^=(a s +c s )-(a s --c s ) cos {9+9') . . 
2<=(a s +c s )-(a s -c s ) cos (O-ff). . . 

We require, therefore, the values of 9 and 9\ 

Now 

cos 9= cos OL cos LN-— sin OL sin LN cos x 

= cos OL(cos LN— - tan OL sin LN cos x) 
cos OL 



« • 



• (i) 

• (2) 



cos X 



cos (LN+X) 



• * 



(3) 



if 



Similarly 



tan A.= tan OL cos x 



si. COS \J xj /T - -.-r /x 

cos 9 = — cos (LN— X ) . 

cos A v 



. . . . (4) 



where 



tan X'= tan O'L cos x- 



Thus to use the formulae for v l9 v % we require to know accurately the position of the 
optic axes. 
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These can only be determined by finding the values of a, b, c. 

The results of the experiments enable us to do this with considerable exactness. 

Let fjb ai fjb h , fju c be the principal refractive indices, /x^ being the greatest, /x^ the least 



P 



a' 



pb 



fr 



1 

C 

1 

b 

1 

a 



Section V. — Values of the Principal Refractive Indices. 

We are now in a position to determine from our experimental results the values of 
a 9 6, c the principal velocities. 

The plane P B,' cuts the plane A C at a short distance from C. 

The section of the surface of wave slowness by the plane A C consists of a circle 
of radius fi b and an oval curve of axes fi a and [i c in the directions OC 3 OA respectively, 
which for a small distance on either side of a principal axis may be treated as an 
ellipse; as in Fkesnei/s theory. 

Also, if L be the point in which P B/ cuts the plane A C 

LP=35°0' 19" 
[vide Section II. (16)]. 

The value of p corresponding to a value of <£'=35° 0' 19" in Table I. for the inner 
sheet will thus be the value of ^ 

But we have as corresponding values 

34°28 / 43 // 1-68115 



34° 53' hT 1-68114 

35° 16' 20" T68115 



we may therefore take 



^=1-68115, 



The values of the radius vector to the outer sheet in the same neighbourhood are, 
Table I, lines 4 and 5, 

34° 47 ; 54" 1-68546 

35° 10' 1" 1-68553, 

For 

f =35° 0' 19* 
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we may take therefore 

/x,= l-68550 

This then must be the length of the radius vector of an ellipse, axes \L a and /x c> which 
passes through the point L — i.e., which is inclined at an angle CL=1° 21' 42" to the 
axis [jL a [vide Section II. (15)]. 

Now if 6 be the angle which a radius vector r makes with the axis /x« 

1 ___ cos 3 sin 3 
r 3 ~~ ^ ft;? 

in the case considered 

0=1° 21 ' 42" 

sin 2 0= '000565 

and the error made in giving to /x c an approximate value will be inappreciable. 

Let us take, therefore, for \l c the value found in the previous part which agrees with 
Rudberg's value 

^=1-53013 

sin 2 



a 2 



= •000241 



Now 



r=l-68550 

-=•351919 

I 

1 / 1 sin 2 0' 



whence 



/jL a ~ cos 2 0\r 2 & 



^=1-68560 



It now remains to determine /x c accurately. 

The natural faces a / m/ (fig. 1) formed a prism whose edge was parallel to the axis O C. 

By passing light through these faces and observing the incidence and deviation as 
before, I was enabled to calculate /x c . 

To pass light through these faces it was necessary, however, to polish them, and this 
operation altered their position appreciably. The values obtained by direct measure- 
ment required a small correction, and to calculate this it was necessary to determine 
the position of these new faces. 

Let the faces approximately coincident with a / and mf be called X Y respectively. 
On levelling the faces of reference m m', I found that Y remained very nearly in the 
same zone as before, viz. : that of m m m / m/ ; while X fell rather below this zone. 
Observations on the angles between X and Y and the faces m m respectively gave 
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XM(= 

XM'(: 

arc YA'M (: 

arc YAM'(: 



= 6> 4 say) 
- B\, say) 
--0 5 say) 
= 6 say); 

Kg. 8. 



123° 32' 
59° 52' 

117° ir 

181° 3' 




Since Y, M, M', are in one zone with A and B, and B bisects arc M M' 

2BA'Y=M'AT+MAY 

.-. arcBAY=]49° 10' 
A'Y=149° 10' — 90° 
= 59° 10' 

Let a 4 /J 4 y 4 be the direction angles of X. 
Then, as before, [Section II. (l)] 



cos a 



v 



sm -^-z — - sin * 



COS fj, 



&c. 



whence 



a 4 = 3° 28' 50 
/3 4 =91° 42' 



// 



To test these values let us calculate the angle X Y ; the triangle X A B gives 



whence 



sin AX 



tt-XAB=60° 43' 50" 



In the triangle X A Y we now know A X, A Y and the angle X A Y. 
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Solving, we get 

XY=119° 5' . . . . (1) 

Experiment gave 

119° 4' 

Here again we have considerable confirmation. 

Fig. 9. 




Let the plane Y X L be the plane of the prism cutting C A C in L ; it will be 
necessary to find A L, X L, and the angle A L Y. 
The right angled triangle LAY gives 

sin AY= tan AL cot AYL. 

The triangle XAY gives 

. T . v a __ s i- n AX sin YAX 

SID. XJ JL xjl ~ ^ rTT 

smXY 

whence 

LYA=3° 28' 30" 

Substituting this value, we get 

AL=2° 59' 10" . . (2) 

From the same triangle we get 

cos YL = cos AY cos AL 

whence 

YL=120° 47' 10" 

but 

XY=119° 4' 

From the same triangle again 

cos ALY= tan AL cot LY 

whence 

F ALY=91° 47" 

2 Y 2 
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Now we know that if / 0' are the points where the optic axes cut the unit sphere 
K, the point where any wave normal meets the sphere if 

KO=6 

KO'=0' 

.•.-^(l+cos^-f?)) 
/v 

=^-j-^-cos(0-6')) 

^ C0B { 2 )~^~^ sm \ 2 ) • •••••• w 

If then we can determine 6— ff since jjl is given by experiment, this formula deter- 
mines the value of [x c . 

Now along the plane X Y, 6— 0' is very small, and a considerable change in the 
positions of O' will produce but small variations in the value of 6—& and .'.of /v 

Let us, therefore, take the value of the angle between the optic axes given by the 
values already determined for fi a and jx h , viz. : 

^=1-68560 
^=1*68115 

and the approximate value of (jl g 

^=1*53013 

If 2<f> is the angle between the optic axes 



tan <A'==— a / J 






whence 

f =9° 4' 5" 

Hence 

AO=AO'=80°55'55" 

AL= 2° 59' 10" 

OL=83° 55' 5" 

0'L=77° 56' 45" 

X=ALK = 88° 13' 

XL=1°43 / 10" 
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Let K be on the side of L remote from X. Experiment shows that this is the 
actual position of the wave normal in the case considered. 
Let 

KX=i// 



.-. LK=KX-XL=f -1° 43' 10" 
cos #= cos OL cos LK+ sin OL sin LK cos ^= -— cos (LK— X) 



cos A, 



tan X= tan OL cos x 
Substituting the values of O L and x 



o i 1-7/ i rjr 



X=16° 17' 10 



cos OL — _ 
log —=1-0427912 



LK-\=t//-18 0' 20" 



Similarly 



PO^J III 

cos 0'= — cos (LK+X') 



cosX 



where 



tan X'= tan O'L cos x 
:. X'=8° l7 / 30 ,/ 

, cos O'L - 

lo g rr= 1-3243204 

° cos X 

LK+X'=f +6° 34' 20" 



A table on the next page gives the values of t/r, the angle made by the incident 
wave normal with the normal to X; D+i the deviation + the angle X Y; i// the 

angle made by the wave normal in the crystal with the normal to X ; and — p being 
the reciprocal of the wave velocity. 
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Table to determine the Value of fi c . 



+• 


T> + i. 


V- 


l 

« 


O / // 


O 1 II 


o / // 




51 16 7 


101 47 15 


30 39 


•427076 


49 41 7 


101 48 45 


29 53 13 


•427061 


47 49 7 


101 57 45 


28 57 43 


•427049 


46 5 22 


102 13 


28 6 3 


•427039 


44 26* 7 


102 34 45 


27 14 5 


•427029 


42 55 27 


103 3 5 


26 25 31 


•427018 



Substituting the values of i// from this table in the formulae for cos 6 and cos & we 
can find 6 and 6\ 

The values oiO—6' are 



o 

3 


31 


50 


3 


37 


30 


3 


42 


20 


3 


47 


50 


3 


52 


50 


3 


57 


40 



Substituting these values in the formula for ^ c , the mean of the results agrees 

closely with 

^= 1-53013 

Thus we may take 

/t„= l-68560 

/ai= l-68115 

^=1 1-53013 



And the angle between O C and the optic axes is 



9° 4' 5" 



e # ■» 



3 <•**«•• 



(5) 



If </> is the angle as seen in air 



O o-i ' er n" 



^=15° 2r 50 



The value Kibchhoff found by experiment was 

15° 27' 

The close agreement is noteworthy, 

As in the first part, the agreement is much closer than it would be taking Rudberg's 
*alues of (A a , fjifo [i c : 
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Section VI. — Application of the Theory to Experiment. — Tables giving Theoretical 

Values of ji. — Discussion of Results. 

We can now continue the work of calculating the theoretical values of the reciprocal 
of the velocity. 

Consider first the plane P B/. 

With the same notation we have [Section IV. (3 and 4)]. 

Fig. 10. 



and 



Whence 



■ 



cos 6= cos ■ cos (NL+L) (1) 

cos A, v y v ' 

tan X= tan GL cos ALN 
cos ff= ^^ cos(NL-X') (2) 

COS A> 

tan X'= tan O'L cos ALN 

OC=0'C=9° 4' 5" [Section V. (5)] 

LC=1° 21' 42" [Section II. (15)] 

.\OL=10 o 25' 47" (3) 

01=7° 42' 23" (4) 

ALN= X = 59 ° 20 ' n " [Section II. (17)]. 

X=5°21' 47". . .' (5) 

log ?^=l-9946698 (6) 

° cos X v ! 

* 
A'=3° 56' 52" (7) 

log— -^= 1-997091 5 (8) 

° cos X v } 
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LN=:LP-PN=LP~f=35°0 / 19"-f (9) 

LN+X=40°22'6"-f (10) 

LN-X'^31 3' 27 // -f (11) 

Substituting these values we can calculate 9 and 9\ and hence find the values of 
v 1 and v 2 or their reciprocals /^ and /x 2 . 

Tables V. and VI. give the results of this substitution : on the right hand for the 
inner sheet, on the left hand for the outer. 

The centre columns give the values of L N+X, L N— X\ 

Then on either side, 6-\-ff for the inner sheet, and d—ft for the outer. Then the 
theoretical values of /^ /x 3 , then the experimental values taken from Tables I. and II., 
and finally, in the outside column, the excess of theory above experiment. 

The next step is to discuss the theory for the plane P Q. 

With the same figure and notation as before, if this plane cut the plane A C in 1/ 
we have [Section IV. (3 and 4)] \ / \' \ being values of X X' in this case. 

Fig. 11. 




whe 



re 



whence 



cos0=^^'cos(L'N+X,) (12) 

cos \ / x " v ' 

tan \ = tan 01/ cos AI/N 
cos ff= co ^ / C os (L'N- X/) (13) 

cos \, x ' v ' 

tan X/= tan O'L' cos AL'N 

00=00=9° 4' 5" [Section V. (5)] 

CL'=1° 15' 40" [Section III. (7)] 

.-. OL'=10°19'45" (14) 

0'L'=7° 48' 25" (15) 

AL'N=x'=59° 13' 2" [Section III. (8)] 
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X /= =5° 19' 45" .......... (16) 

X/=4° 0' 50". . (1^) 

L'N=L'P+PN (N being now on the side of P remote from I/) =35° 3' 14"+<£'. 

L'N+X / =40° 22' 59"+(// (18) 

I/N-\/=31° 2' 24"+f (19) 

l og 5^2i [ i = ^9947 864> ........ (20) 

& cos\ v ' 

r0 Q O'T/ - 

log __ —i-9970226. ........ (21) 

and as before, we can calculate from these data the values of 0, & and hence of ft x /%. 

Tables VII. and VIII give the results of this. 

The fifth column gives the values of L N+X for the inner sheet, the next L N — \' 9 
the next 0+0', the eighth the value of the reciprocal of the velocity calculated on 
Fresnei/s theory [by means of Section IV. (1 and 2)], 

— i=a *+c*-(a*—c?) cos (0+0') 
Pi 

^ =a »+ c 3_( a a_ c a) cos (0-0') 

the ninth, the same quantity from experiment, and the last, the excess of the 
theoretical over the experimental value. 

The values of 0--0 / given in the fourth column are calculated for the inner sheet. 

Now as the values of <{/ corresponding to the same incident wave are different for 
the two sheets, the values of and 0' will be different ; but on referring to the values 
it will be seen that after the first twenty observations, 0— 0' is very nearly constant 
throughout the arc considered ; so that /x 2 varies very slowly, and except just at first, 
we may, to the degree of approximation required, use the value of and & for the 
inner sheet in calculating the theoretical value of the reciprocal of the velocity for the 

outer sheet. 

In the case of the first twenty observations for which 0-— 0' varies appreciably with 
<£', the experimental values of /^ are obtained by interpolation from those values in 
Table I., to which they correspond. 

Table V. contains the theory for I. 

„ vii. „ „ in- 

„ VIII. „ „ IV. 
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Let us consider first the inner sheet. 

In Table V. we see the differences are small to begin with, and increase fairly 
regularly throughout the whole of the table, the experimental value being almost 
always less than the theoretical. 

Lines 4, 5, and 6 were those chosen to determine the value of /x# which explains the 
coincidence of the two curves in that neighbourhood. The table covers an arc of the 
surface of nearly 20°, as is seen by referring to the values of <fi f in Table I. 

Throughout this arc the theoretical section lies outside the experimental, the 
difference between the two increasing as we proceed from the axis //,#. 
This continues throughout Table VI. 

The first three lines of Table VI. just overlap the last three of Table V. The difference 
reaches a, maximum value towards the end of Table VI. at a point for which <fi f = 
1° 19' 54" (vide Tables II., 21, and VI., 21), that is at about 34° away from the axis /x^. 
After this the difference begins to decrease at first slowly throughout Table VII. 
which refers to the other prism P Q, so that the arc of section has changed slightly 
between Tables VI. and VII., both arcs passing through the point P. 
Table VII. covers an arc of 22° (Table III.). 

The first two lines in Table VII. refer to a wave normal lying on the same side of 
P as L. In the rest of the table the wave normal is on the side of P remote from L. 

This arc is continued throughout Table VIII. ; the theoretical values of /x are still 
greater than the experimental, but the difference diminishes as we proceed along 
the arc. 

The last three observations in Table VII. and the first three in Table VIII. overlap. 
Table VIII. covers an arc of 20°. 

Thus the observations extend over two arcs, each passing through the point P 
inclined to one another at that point at so small an angle as to be almost continuous, 
the one arc being 35° 40', and the other 40° 40' in length. 

The results of theory and experiment agree at the extremity of the first arc, differ 
most widely in the neighbourhood of the point where the two arcs meet, and tend 
towards equality again throughout the second arc. 

It is worth noticing that the experimental results at the end of Table VI., the end 
of the first arc, agree closely with, those at the beginning of Table VII,, the beginning 
of the second arc. 

The experimental results for the inner sheet are therefore represented by a curve 
which coincides with Fresnel's curve at the extremity of one axis, and lies inside it 
throughout the next 75° of its length. 

For the outer sheet the differences are less than for the inner sheet. 
Q—& varies so slowly that it seemed sufficient for the purposes of comparison to 
calculate every fifth theoretical value, with the exception of that portion of Table V. 
for which the variations of 6— 6 r are sensible. 

The section cuts the principal section A O C between lines 4 and 5, Table V. Just 
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at that point the theoretical value is greater than the experimental. It continues so 
until we reach line 18, Table V., that is for an arc of about 9°. 

From that point onwards the theoretical value of /x 3 is uniformly less than the 
experimental. 

Thus the result of experiment for the outer sheet is represented by a curve which 
coincides with the theoretical section at the point in which it cuts the plane A O C, 
then lies inside it for an arc of about 9°. 

At this point the two again coincide, and for the rest of the arc experimented on, 
the section given by experiment lies outside that given by theory ; the difference 
between the two increasing throughout Tables VI. and VII., and diminishing again in 
Table VIII. 



Section VII. — Effect of possible Errors in the Values of \L a > /x s , ji c discussed. 

The next task will be to determine the effect of any small errors which may have 
occurred : 

(1.) In the determination of \x a , /x#, \x c \ 

(2.) In the angles which fix the position of the planes P R, P Q relatively to the 
axes of the surface. 

(1.) fM a . To determine \i a we considered the length of the radius vector common to 
the principal section A O C and the section by the plane P R. This radius vector is 
fixed in length, and is on Fresnei/s theory the radius vector of an ellipse, axes \i a and 
He inclined, at a small angle to /x^. 

The increase in the value of \x a will therefore increase the value of this angle, that 
is, it will alter the position of the plane P L R' and be more properly considered 
under (2). 

\li is determined by the intersection of the section of the inner sheet by the plane 
P R, with a circle of radius \x h and is therefore fixed at least within the limits of 
experimental error. 

fM c is found from an independent and, on the whole, less trustworthy series of 
observations — less trustworthy because the faces were much less plane than P, Q, R. 

Let us consider the effect of decreasing the value of (jl c . 

If e be the angle between the optic axes 



tan e: 




As fx c decreases tan e decreases; therefore € decreases. 
Let us consider how this affects and ff. 
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Taking the same figure as previously. 



Fig. 12. 




Let O' K be an arc perpendicular to A C. 

Let J 0/ be the new positions of 0'. 

Then 

0,0:= 0/0'= Se 

Let 0, m 9 O' n be perpendicular on N 0, N 0/ respectively. 

OX =7° approximately 
O'LP=60° approximately 
.'. LK=14° approxinmtely 



Let N fall between P and K 

80= -Om= -8 € cos NOL 

S0'=O>=SecosNO/A 

8(0+0')= -8e (cos NOL- cos NO/A) 

Thus 8(0+ 0') is negative and decreases numerically as we approach P ; for N L, 
N 0/ A become more nearly equal. 

Thus between K and P 6-\-6' decreases with e. 

Between L and K both and 0' decrease, therefore 0+0' decreases d fortiori; but 
8(0+0') is always numerically less than 2Se. 

Let us see how these changes affect /x x . 

We have [Section IV. (1)] 



2 



1_ _1_ 

2 • TTs 



1 



Ah fta ftc 



,/V* Pa 



cos(0+0 / )« 



11.. 11.. 
If fjb c is decreased, ~^+— is increased, —% ^ is increased, 0+0' is decreased. 

Therefore cos (0+0') is increased so long as 0+0' is less than 90°. 
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Therefore f — A cos (0+0') is increased so long as 0-\-ff is less tlian 90°. 

And the increase in ( — ^— —) cos {0-\-ff) is greatest when 0-\-& is least. 



J^c ^ 



a 



Also the increase in l— -] cos (0+0') is less than the increase in ( -z+ — 

Therefore — -J — - — (— 3 3 )cos(0+# / ) is increased, and this increase is greatest 

f^c fta Xftc fat J 

when 0-\-0* is greatest. 

Therefore, as 0-\-0' increases up to 90° \l x decreases. 
When 0+ff is greater than 90°, cos (0-\~0 f ) is negative. 
As 0-\-0' increases, cos (#+#') increases numerically. 

Therefore ( — A cos (0+0') increases numerically and is negative. 

1 1 / 1 1 \ 

Therefore — ^\ — 3 — \~ % —~~A cos (0+0') increases a fortiori. 

Therefore \l y decreases. 

Thus the effect of decreasing \l g is to produce a decrease in the value of [i 19 and this 
decrease continually increases as we go from L to Q. 

Now if we compare the differences between the theoretical and experimental values 
of [a 1} given in Tables V. and VI., we see that they increase with 0+0', the theoretical 
value being in excess of the experimental. 

The effect of the assumed alteration in the value of \i c would be to decrease the 
theoretical value of /x x by an amount which constantly increases, that is, to bring it 
more nearly into agreement with experiment. 

But the effect in Tables VII. and VIII. would be contrary. 

On referring to them we see that the differences between theory and experiment 
continually decrease. 

The effect of the proposed alteration would be to subtract from the theoretical value 
of /x x given in those tables : a quantity which increases rapidly as we get further 

from L. 

The result would be that the theoretical value would soon become less than the 
experimental, and the difference between the two would continually increase. 

We might, it is true, choose the decrement of \i c so as to make the extremities of the 
arcs considered coincide, and then the differences would be diminished throughout. 
The effect of this would be to make the theoretical section lie outside of the experi- 
mental for about 75°, then cut it and he inside it. 

Whether the two could again be made to coincide when cutting the plane A O B 
experiment alone could settle. 

Let us now consider the effect of the supposed change on /x 3 . 

We have 

MDCCCLXXIX. 3 A 
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2 ^-jWA-AW*-*'). 



^2 f^a l^G \f^c f^a 

We have seen that is decreased by the change. 

& is at first decreased, though by less than 0, afterwards it is increased. 
Therefore, for any point on the curve L P the value of 0—0' is decreased by this 
alteration. 

Therefore cos (0-—0 f ) is increased. 

( — -) cos (0—0') is increased, though by less than ( ~^+~^). 

1 1 / 1 1 \ 

Therefore ~^+ - ~ 3 — " 3 {] cos (0—0') is increased. 

pc pa \pc pa J 

Therefore \l % is decreased. 

Or the variation in p% is in the same direction as in p v 

But since 0—ff is <6+ff (~*~~*) C0S < < 6 ~ ff ) is > fV~A) cos {0+&)- 

And the whole change in 



1/1 1 



pc pa \pc pa 



cos {d-ff) 



is much less than the whole change in 



\+-^-l 1 -- 9 )cos(0+0y 

spc" Pa) \Po pa 



Or the decrease in /x 3 is small compared with that in [i^ 

Again, as the variation in 0—0' varies from zero at L to a maximum about K, while 
cos (0—0') also varies considerably, the variation in /jl 2 will be considerably different 
for different points between L and K. 

After passing K, however, 0— & is nearly constant all along the curve. 

The change in the value of fi 2 , therefore, will not differ much for different points, 
that is, for different waves. 

The whole effect of the alteration, therefore, will be to decrease throughout the 
theoretical value of /x 2 , this decrease being greatest about Table V., line 23. 

In the greater part of Table V. the result will be to produce a closer agreement 
between theory and experiment. 

But in the rest of the work, since the theoretical value of /x 2 is less than the experi- 
mental, the effect will be to increase the difference between the two and so widen the 
discrepancy between theory and experiment. 

Thus taking both sheets we cannot produce on the whole a closer agreement 
between theory and experiment by decreasing the value of [i c . Neither can we by 
increasing it. 

For though we might produce closer agreement for the outer sheet by thus increasing 
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the theoretical value of ^ s , we should acid greatly to the error in the case of the inner 
sheet. 



Section VIII. — Effect of Change in Position of Planes of Prisms with reference to the 

Principal Planes of the Crystal. 

We have now to discuss the effect of small variations in the position of the planes 
P II and P Q. 

We take P R first, and we consider its position as determined by C L and the 
angle A L P, while the position of the normals to P and R is fixed by P L (vide 
Section II.). 

Eig. 13. 




Now variations in C L and ALP produce variations in L, 0' L, X X', &c, and 
hence in 6 ff 9 while variations in L P only change 6 and ff, L, O' L, X and X 7 remaining 
constant. 

Let us consider, therefore, the variations of L P first. 

A reference to our tables shows us that through the greater part of the arcs L P, 
P Q', 0—ff is constant, so that §6 considered as dependent on L P is equal to SO' 
nearly. 

Again, 



Pi \pa /VV \M 



1 +A N )-( i 5 - i - a )cos^-^) 



Ms ^M« M<? / Xftc fta 



and since the change in (6—9') due to a small increment of L P is very small, such an 
increment will not alter the value of p 2 . 

But if 6—ff is nearly constant for a small motion along L P, 6+ff will be nearly so 
for a small motion perpendicular to L P. 

Thus a small motion perpendicular to P L will not alter appreciably the value of ^ ; 

o .A. j-i 
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we may therefore obtain corrections to [i x /x 2 independently, to /x x by supposing an 
increase in the value of L P, to /% by supposing the planes P R. P Q to turn through 
a small angle about some line in them. 

Moreover, we shall consider that line to be in each case the line of intersection with 
the plane AOC, and treat this line as the same for both planes, the two lines thus- 
considered as coincident are really inclined at an angle of 3'. 

To find then the variation which we must assume in L P to bring the results of 
theory and experiment into closer agreement for the inner sheet, we proceed as follows — 

If in the equation 

A = /± 3+ J_\_(i__-l\ cos (Q+ff) 

we substitute an experimental value for p l9 we can determine a value of {0-\-0'). 
Subtracting this from the calculated value, we get B(0+6 / ). 
But for the displacement considered 80=80'. 

Hence we have found 80 : but cos (LN — X) = -p-=r- cos 0. L and X are unaltered 

v J cos OL 

by the displacement supposed; therefore we get 8(L N), and LP=LN+f. j> is 
constant, therefore we have found 8 L P. 

To put this plan into execution, I chose as the experimental value for ^ the 21st 
in Table VI. 

/ x 1 =l*62807 

the differences between theory and experiment having at that point a maximum value 
and being fairly regular. 

The resulting value for (#+0 / ) is 

70° 57' 20" 

The calculated value is 

70° 16' 40" 

The difference is 

0° 40' 40 

.-. 80=0° 20' 20 
This gives for S(L P) 

8(LP) = 0° 21' 20" 

Thus at this point the result of theory and experiment could be brought into 
coincidence by supposing the line O P normal to the face P to turn through an angle of 
0° 21' 20" in the plane L P away from L. 

Before discussing the effect of this change on the values of /x x in other directions, 
let us consider the variation it will be necessary to make in the angle A L N to bring 
the theoretical results for the outer sheet into accordance with experiment. Taking 



// 



// 
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the experimental value for /x 2 corresponding to the incident wave which gave rise to 
the value of /x x we have just been discussing, we have 



^=1-68447 



Whence by means of the formula 

we find 

0-0'=9° 15' 

The value of 0—0' given by theory is Table VI., line 21, 

6- 6' =9° 30' 

We have to find the change in the value of x or A L N which will produce this. 
In this variation 0-f-0' is constant 

.-.80= -80' 
but 

8(0-0') = -15' 

.-. 80= -7' 30 
86'= Y 30 



// 



// 



From these results we can find the increment in x from the equations 

cos 6= —-cos (LN— X) 

cos X N 

tan X= tan OL cos x 

cos 6=. cos OL cos LN+ sin OL sin LN cos x 

^ sin 080 

^ sin OL sin LN" sin % 

> Substituting values we are led to a value of S x> near ty equal to 1°. So that the new 
value of x =60° 20' [Section II. (17)]. 

This change in the value of x produces a corresponding change in the values of X, X'. 

The new values are for the prism P II 

X'=3° 50' 0" (2) 
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while 

log G °^= 1-9945594 ......... (3) 

COS A 



for the prism P Q 



cos OL 
eosV 



1-9970325 ......... (4) 



A — JL U Z U . . . « . . , . . * . (/ 

A '— O 41: o . . , , . , . . , ,, . \v) 

log CQS =1-9946764 ......... (7) 

° cos X / v 7 

io2f r j — — 1 yyoyozy . <> . » , , . . , (q) 

° cos A, v ~ ! 

The ensuing Table gives the values of fi l9 /jl 2 for nine different incident waves taken 
from the previous work, calculated on the supposition that both these changes have 
been made in the position of the plane P P. 

The column headed Difference gives the excess of theory over experiment ; that 
headed a gives the excess of the theory before it was modified ; /3 gives the change in 
the theoretical value of fi x due to this displacement. 
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It will be noticed that the differences for this modified theory are much less than 
for the original, but they have nearly all changed sign. In fact, the correction applied, 
calculated so as to produce the closest possible agreement in line 5, has been too great, 
and the differences can be reduced by applying a similar but smaller correction in the 
opposite direction. 

Now this variation in the value of /^ has been caused almost entirely by the 

increase of the angle LP; we must therefore decrease this angle again, and as the 

numbers in the difference column are each about one-fifth of the corresponding number 

in the column /3, we must decrease L P by one-jifth of the amount by which it was 

previously increased. 

But this increase was 

S(LP) = 0° 21' 20" 

we must therefore decrease this by 



leaving the total increase 



0° 4' 15' 



0° 17' 5' 



In consequence of this, we must add to the modified theoretical values of /^ 
respectively, quantities equal to one-fifth of the corresponding numbers in column /3. 

The result gives us a new set of theoretical values, which agree with experiment 
throughout to a remarkable degree. 

These results are expressed in Table X, 

Column 1 gives the values of ^ already modified. 

Column 2 gives one-fifth of quantities in /3, which, when added to the values of 
/x x in Column 1 respectively, give Column 3 the resulting theoretical values. 

Column 4 gives the experimental values, and Column 5 the excess of theory over 
experiment. 

""able 



Modified value 


~5 r * 


Final value of /x^ 


Experimental 
value of ^ j. 


Difference. 




1-68102 


1 


1-68103 


1-68099 


4 


1 


1*67707 


7 


1-67714 


1-67721 


-7 


2 


1-66285 


13 


1-66298 


1-66300 


-2 


3 


1-64590 


17 


1-64607 


1-64603 


4 


4 


1-62806 


18 


1-62824 


1-62807 


17 


5 


1-60882 


18 


1-60900 


1-60897 


3 


6 


1-58346 


17 


1-58363 


1-58365 


—2 


7 


1-55138 


16 


1-55154 


1-58157 


-3 


8 


1-53779 


5 


1-53784 


1-53774 


10 


9 



These differences are exceedingly small— sometimes positive, sometimes negative- 
but rarely greater than the possible error of experiment as appears from the tables 
of error in I., II., III., and IV. 
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Thus we have found a plane section of the surface which agrees closely with the 
results of experiment. 

This section passes through L, a point on the principal section A C at an angular 
distance of 1° 2 1 7 42" from C, is inclined to that section at an angle of 60° 20', while 
the normal to the face P lies at an angular distance from L equal to 



35° 0' l^+O 17' 15"=35 6 17' 34 



// 



Under these circumstances it will be necessary to consider the probability of making 
such an error as this result indicates in the position of the normal to P with reference 
to the principal axes. 

But before doing this we must say a few words about the variations in the values 
of /x 2 produced by the correction now considered. 

The final correction will not affect any of them except the first, which will be slightly 
increased by it. 

This follows from the fact that between lines 1 and 2 of Table IX. 6 — Q' varies 
considerably, and increases or decreases with #+#', so that fi l9 /x 2 increase or decrease 
together. 

The change considered has increased /x D it will therefore increase /x 2 * 

The difference in line 2 is considerable, but it may be noticed that the experimental 
value is clearly too small, falling as it does (Table I., line 18) between two con- 
siderably greater values. In the other cases the differences between theory and 
experiment have been diminished but remain of the same sign as before. 

A still further displacement in the same direction would therefore produce still 
further agreement. 

We have now to consider the effect of Varying C L< 

Let us trace the effect of increasing C L by 10'. 

The value of fjb a depends on that of C L being found from it by the formula 
(Section V.) 



--== sec^UL 



/V \1 S P? 



where 



p*= 1-53013 



and r is given from experiment 



r==l-68550 

The new value of a will be found to be 

/v=l-6857i 
instead of 1*68560 

This increase in the value of a will increase the angle betwen the optic axes from 
9° 4' 5", the value already found, to 9° 9' 20". 

MDCCCLXXIX. 3 B 
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We have, therefore, to discuss the effect of these changes in the theoretical values 
of p x , /x s . 

We will further suppose that the angle \ is so altered that the section still* passes 
through the former position of P — i.e., so that the only error made in the previous 
determination of its position has been in the position of R'. 

This is probable, for C B,' is small and A R/, B IT, nearly right angles. 

The new value of x W U1 b e 

x =59 6 32' 

* 

These changes will produce variations in the values of L, (XL, X X', for P R/, and 
O I/, 0' I/, X / X/, for P Q (we assume that L L' is unaltered). 

The complexity of the changes rendered it very difficult to discuss the effects on 
ju, 1? fa from general considerations. I have therefore calculated the new values in five 
different positions taken from the previous tables* 

The table following gives the general results, 

Table XI. 



Whence taken. 


Difference. 


Mi 

Modified Theory. 


Mi 
Experiment, 


Mi 

Former Theory* 


Difference. 


1 

2 
3 

4 
5 


Table V., line 28 
VI 11 

,5 * x,, ,, XX 

„ VI., „ 21 

* > T Xl.» A * M _X. 

^ * T JL JL JL a a a a _X. -L. 


36 
49 
65 
41 
43 


1*66386 
1-64652 

1-62872 
1*58406 
1-55200 


1-66300 
1-64603 
1-62807 
1*58365 
1-55157 


1-66352 
1-64674 

1-62897 
1-58432 
1*55216 


52 
71 

90 
67 
59 



Whence taken. 


Difference. 


M2 

Modified Theory* 


M2 

Experiment. 


M2 

Former Theory. 


Difference. 


1 

2 
3 
4 
5 


Table V., line 28 
„ VI., „ 11 
„ VI., „ 21 
„ VII., „ 21 

ft a Y JL JL JL. • a a a _L _L 


-1 

-6 

+ 1 

-7 



1-68447 
1-68446 
1-68448 
1-68450 
1*68452 


1-68448 
1-68452 
1-68447 
1-68457 
1-68452 


1-68437 
1-68437 
1*68437 
1-68439 
1-68440 


-11 

-15 

-10 
-18 
-12 



The difference table on the left is between the modified theory and experiment, that 
on the right between the former theory and experiment. 

The upper table refers to the inner, the lower to the outer sheet. 

It will be seen that the differences for j-x x> /x 3 are both decreased, though still 
remaining of the same sign. 

The differences for /x x are about two -thirds of their former value. 

If we then increase C L still further by twice as much as previously, so that the 
total increase is 30', we shall get still closer agreement for /x ]B 



MR. R. T. GLAZBBROOK ON PLANE WAVES IN A BIAXAL CRYSTAL, 



371 



Table XII. 



a. 


j8. 


Final Theory. 


fx 2 
Experiment. 


Difference. 


16 


48 


1-66304 


1-66300 


4 


22 


66 


1-64608 


1-64603 


5 


25 


75 


1-62822 


1-62807 


15 


26 


78 


1-58354 


1-58365 


-11 


16 


48 


1-55268 


1-55157 


11 



r- 


d. 


t*2 

Final Theory. 


A*2 

Experiment. 


Difference. 


10 


30 


1-68467 


1-68448 


19 


9 


. 27 


1-68464 


168452 


12 


11 


33 


1-68470 


1-68447 


23 


11 


33 


1-68472 


1-68457 


15 


12 


36 


1-68475 


1-68452 


23 



Column a gives the change in /jl 1 due to a change of 10' in C L. 

Column /3 gives three times that change which, assuming the theory of proportional 
parts, will be the change in /x x due to a change of 30' in C L. 

Subtracting this from the values of fi Y given by the original theory, we get the final 
theoretical values of fi v The$e agree much more closely with experiment than the 
original values. 

But the effect of this further change will be to over correct /x 3 . 

The first change has reduced the differences for /x 2 by the amounts given by y, 
Table XII. 

Multiplying these by 3 we get S, the changes produced in /jl 2 by a change of 30' 
in C L, and adding these values to the original theoretical values of /x 2 we get the 
final theory, which differs from experiment by rather more than the original theory 
and in the other direction. 

Thus an increase of C L will not, on the whole, produce the required effect. 
Similarly a decrease will not produce it either. 



Section IX..— Possibility of an Error in the Positions of the Faces discussed. 



We have seen, however, that by increasing LP by 17' 5" and A L P by 1°, we can 
make the agreement between theory and experiment extremely close. 

This alteration in the values of L P and ALP may be effected in two ways. 
We may suppose either that, the axes of elasticity retaining their position with 

3 B 2 
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reference to the planes m / m\ we have made an error in the observed angles P M, P M' 
(Section II.), and hence in the position of P referred to the axes. 

Or, secondly, that the position of P relatively to m m has been accurately found, 
but that the axes of elasticity differ slightly in direction from their assumed positions, 
so that O C is not parallel to the intersection of m m, and B does not bisect the 
angle M M'. 

We will take the first hypothesis first, and consider what must be the error in P M, 
P M' to give rise to the required change in the position of P. 

Fig, 14 




Let P' be the new position of P, we require M P' M' P\ 
Let P K be perpendicular on L P ' 

LK=LP approximately 
KP'==0° 17' approximately 
KP = PLK sin 35° approximately 

= 60 / X*57 

= 34' approximately 

Now P K P / is approximately a small plane triangle with a right angle at K. 



PP'^VPE^+P'K 3 



whence 



Again 



whence 



V / 34 S +17 3 



X X — o o 



cos CP'= cos CL cos LP'— sin CL sin LP' cos ALP' 



CP'=35° 58' 

CP =35° 43' [Section II. (12)] 

PF= 0° 38' 
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Hence from the triangle POP' we get 

PCP'=Q° 58 
But from the triangle A C P 

ACP=57 42' 30" [Section II,] 

s \ ACP'=58°-40' 30" 

but 

ACM=58° r 30" [Section II. (7)] 

.•. MCP= 0° 33' 

cos MF= cos MCP sin CF 

whence 

MF==54° 58' 50" 

The original value of M P / was as given by the mean of a large number of 

observations 

MP'=54 G 17' 6" [Section II, (3)"]. 

The difference, amounting as it does to over 40', is far beyond any possible error of 
experiment. 

Let us find further the change in M' P, 

cos M'F= sin CF cos M'CF 

M'CP'=63° 47'-33' 

= 63° 14' 
whence 

M'F=74°-39' 40" 
The observed value was 

M'P=75° 16' 20" [Section III (4)]. 

The difference is 0° 36' 40". This, again, is far beyond any possible experimental 
error. Thus, so long as we assume the position of the axes of elasticity to be 
definitely fixed in the crystal so that O C is parallel to the intersection of m m while O B 
bisects the angle between them, the displacement of the plane P B/ necessary to bring 
the results of experiment and theory into agreement is far too great to be possible. 

But there remains the supposition that the axes of elasticity have not exactly the 
same position in all crystals of arragonite. So that the displacement of P B relatively 
to the axes might be effected by changing slightly their position with reference to the 
faces of the crystal, P B retaining their position, relatively to those faces, unchanged. 
The possibility of this is a question for the mineralogist. I have been as yet unable to 
find data for a satisfactory answer. It seems, however, plausible to suppose that in a 
substance like arragonite which is not chemically pure, but contains foreign substances 
to a variable degree in different specimens, some slight variation such as that indicated 
might occur. 
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We can, without much difficulty, calculate the amount of the change. 

Fig. 15. 




We may very approximately treat the rotation round O L as if it were about C, 
C L being very small. 

The effect of a rotation of 1° about O C will be to bring A and B to A, B 7 respec- 
tively in the same plane, where 

AA,=BB,=1 

A / B 7 being to the right of A, and B. 

The increment of 17' in L P may be effected by turning the axes through an angle 
of 17' about a normal to the plane L P passing through O, let S T be this normal, 
let A' B' C be the final positions of A B C, 

A, A', B, B' are arcs of small circles whose centres lie on S T at T and S respec- 
tively, so that 

A / TA , =B / SB , = 17 / 

C (7 is an arc of a great circle, and is =17'. 



B/B^SB^B'SB, 

806,= 60° approximately 



SB'= sin 60°= 

.-. b;b=^^ 



o v'3 



BB'=V^£/+B'B/) 

<3xl7 2 \ 



=v / ( 6 ° /3 + 



4 } 

■ 1° 2' approximately 
A.A^ATA'sinSO 

= 8' 30" 



AA' 3 = \/ ( 60' 2 + "77 ) approximately 



1° 1' approximately. 
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Thus the new positions of the axes are inclined to the old at angles of 1° 1', 1° 2', 
and 17 '' respectively for A, OB, 0. So that the axis of in its new position 
would be inclined to the line of junction of m m' at an angle of 17 '', while that of b 
would be inclined to the bisector of these planes at an angle of 1° 2'* 

If variations such as here considered be possible, so that we cannot be certain 
d priori of the position of the axes of elasticity, the only method of testing Fresnel's 
or any other theory will be by triaL 

Section X. — General Results of the Investigation. 

Combining the results of the two series of experiments, it seems to me most probable 
that Fkesnel's theory is only true as a first approximation. 

Both series of observations have led to the discovery of considerable deviation from 
the theory, unless we assume errors in the experiments, especially in the determination 
of the position of the planes of the prisms with reference to the axes of the crystals, 
which are greatly in excess of the amounts we can reasonably expect* 

It will be noticed that the values taken for the constants fi a and //,$ in the two parts 
of the paper differ by *0002 and '00017 respectively. This is due to the fact that the 
two crystals used were different. Rudberg (Pogg. Annalen, vol. xvii., p. 1) found 
differences of as much as '0004 in the values of /x# deduced from two specimens of 
arragonite. 

It may be objected that the variations between theory and experiment are not in 
exactly the same direction in the two Parts ; but we must remember (l) that the 
arcs investigated are taken from entirely' different portions of the surface ; (2) in the 
first Part, the approximately elliptic section belongs to the outer sheet; in the second, 
to the inner sheet of the wave surface* 

Fig, 16, 




The accompanying figure will illustrate the results of the investigation. 
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M N is the principal plane of the first prism nearly coincident with the plane 
through the optic axes A C. 

KOL the principal plane of the second prism nearly coincident with the plane 
B A, so that N is nearly coincident with A. 

M L the principal planes of the two prisms in the second crystal here treated as 
coincident, inclined at about 60° to B C, and cutting it in a line nearly coincident 
with C. 

Hence M is not far removed from C. L is nearly in the plane A B and 
inclined at about 60° to A or O N» 

The strong lines give approximately the form of the sections of Fkesnel's surface 
by these planes, the dotted lines the results of experiment. 

In the case of the arc N' 1/ the results of theory and experiment agreed closely. 

For the arcs M N, M' 1ST the experiments covered an arc of about 16° from M 
and M'. 

For N L the experiments covered an arc extending from L to about 10° on the side 
nearest to N of the point where the two arcs intersect. 

For M L the experiments extended over an arc of about 70° measured from M. 

Section XI. — Effect of Dispersion considered. 

The theory of dispersion appears to me to afford a more probable explanation of 
these small variations from Fkesnel's -construction. 

Fkesnel himself remarked ('Second Supplement au premier Memoire sur la 
double refraction/ GEuvres Completes de Fbesnel, torn, ii.) that in the case of the 
vibrations which constitute light the radius of the sphere of action of the molecular 
forces brought into play by the vibration is not necessarily Very small compared with 
the wave length. 

And, consequently, it is incorrect to suppose that the propagation of each of the 
disturbances of which a vibration is composed is uninfluenced by the disturbances 
which precede and follow it, and that the velocity of propagation is independent of the 
manner in which they proceed and follow it. 

This supposition is the basis of Fresnel's work on double refraction. 

Let us consider the effect of dispersion in a doubly refracting medium. 

In an isotropic medium the relation between V, the velocity of wave propagation, 
and X } the wave length, is generally allowed to be 

a, 6, c, &c, being constant, the values of the terms continually decreasing, so that 

1 b 

except in highly dispersive media we may put - =a+^i with sufficient exactness. 
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Let us suppose that an equation of this form holds in crystalline media also, only 
that a, b, c, &c., instead of being constants are functions of the directions of propagation 
and vibration; and, further, let us suppose that Fresnel's construction is true for 

1 

waves of infinite length, so that the equation -=a gives us a Fkesnel's wave surface. 

From the known values of the constants of the wave surface for different rays of the 
spectrum the constants of the surface for infinite wave lengths can be found, and hence 
the value of a calculated in any given direction. 

1 

If from experiment we find the value of - or /x for any wave length, the difference 

h 
fL—a ought on this theory to be equal to ~z. 

And if we find the values of /x for different wave lengths ((jl 1} /x 2 , /x 3 say) in the same 

direction, we have 

b 



/-&o *"" "" Cl/' 



V 



/x 3 — a — ^ 
whence 

fly — & K 



2 



L6o — ~~ Hi i\){ 



2 



'2 



To verify this I observed the values of /x in two directions for the rays 0, D, and F, 
with the following results : — 






= 1-2403 



fJbjy—CL 



Po 



—a 



/% — a 



1*2875 (first experiment) 
1*2770 (second experiment) 

^Y=l-46978 

1*47208 (first experiment) 
1 1*47348 (second experiment). 



The numbers, especially in the last case, are sufficiently close to make it seem worth 
while continuing the investigations, 

mdcoolxxix. 3 o 



